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We study possible spin liquids on square lattice that respect all lattice symmetries and time- 
reversal symmetry within the framework of Schwinger boson (mean-field) theory. Such spin liquids 
have spin gap and emergent Zn gauge held excitations. We classify them by the projective symmetry 
group method, and hnd six spin liquid states that are potentially relevant to the J 1 -J 2 Heisenberg 
model. The properties of these states are studied under mean-held approximation. Interestingly 
we hnd a spin liquid state that can go through continuous phase transitions to either the Neel 
magnetic order or magnetic orders of the wavevector at Brillouin zone edge center. We also discuss 
the connection between our results and the Abrikosov fermion spin liquids. 
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Quantum spin liquid is the ground state of a quan¬ 
tum spin system in two and higher spatial dimensions 
that does not show any spontaneous symmetry breaking. 
Since this concept was first introduced four decades ago-H 
a lot of theoretical and experimental progress have been 
made in search for this unconventional phase--. In par¬ 
ticular several promising experimental candidates have 
been identified- 2 3 . On the theoretical side, it is un¬ 
derstood that quantum spin liquid is more likely to be 
found in frustrated spin-1/2 models with large classical 
ground state degeneracy, where strong quantum fluctua¬ 
tions within the degenerate classical ground states may 
prevent long-range symmetry breaking ordef^. 

The square lattice Heisenberg model with nearest- 
neighbor ( Ji) and next-nearest-neighbor (J 2 ) Heisenberg 
couplings (abbreviated as J 1 -J 2 model hereafter), 

H=J 1 J2S t -S j + J 2 J2 S ^ S ^ (!) 

(n) am 

is one of the simplest frustrated spin models. It has at¬ 
tracted a lot of attention for its possible relevance to 
the copper oxidetPI 

and iron-based high-temperature 
superconductors(HTSC)P^. Ground state of this model 
is the Neel order [(S^^) oc (—l) a:+y ] for Ji | J 2 1 > 0 
which is the magnetic order of undoped cuprates, or 
the “stripe” collinear order [(S/^.y)) oc (— l) x or (—l) y ] 
for J 2 |Ji| > 0 which is the magnetic order of 
many parent compounds of iron-based HTSCs. Spin liq¬ 
uid physics in this model has also been proposed to be 
relevant to the high temperature superconductivity in 
these materials® 11 . However it should be noted that 
the appropriate model for iron-based HTSC is possibly 
spin-1 (unlike spin-1/2 for cuprates) and may involved 
biquadratic spin interactions 12 . In this paper we will 
mainly focus on the spin-1/2 case. 


* Current address: Department of Physics, Boston College, Chest¬ 
nut Hill, Massachusetts 02467, USA 


The classical ground state of the J\-Ji model at 
J 2 /J 1 = 1/2 has large degeneracy. The ground state 
for quantum spin-1/2 model around J 2 /J 1 = 1/2 has 
no magnetic order, although the true nature of this 
disordered phase has been under debate for a long 
timJHHUl Recently a density matrix renormalization 
group(DMRG) study^ showed evidence of a spin liquid 
phase with spin gap in the region 0.41 < J 2 /J 1 < 0.62. 
However a later DMRG result of the same model did not 
confirm this conclusion®. We do not intend to answer 
whether a spin liquid phase exists in the J 1 -J 2 model. 
The goal of this paper is to classify spin liquid states 
within the Schwinger boson formalism, identify candi¬ 
date states relevant for the J\- J 2 and related models and 
study their properties. It is possible that none of the 
spin liquid states studied in this paper can be realized 
in the Heisenberg models on square lattice. Neverthe¬ 
less they may still be ground state for models close to 
the J 1 -J 2 model, for example models with ring-exchange 
interactions. 

Spin liquids generically have f ractio nalized spinon and 
emergent gauge field excitation d— I— ( The spinons are 
spin-1/2 and can be bosonic or fermionic. In this paper 
we restrict ourselves to the bosonic spinon (Schwinger 
boson) representation. The spin operators on every site 
i are expressed by two bosons as 

Si = \ E b l" a p b W ( 2 ) 

a,0=7,1 

This enlarges the onsite Hilbert space, a local constraint, 

^ = X fe Ua = 25, (3) 

a 

should be enforced and can be implemented by introduc¬ 
ing a U(l) gauge field. In this formalism, spin liquid 
will have gapped bosonic spinon s and boson condensa¬ 
tion will produce magnetic orderdMMI j n order for this 
state to be stable in two spatial dimensions(2D), the U(l) 
gauge field should be gapped, the simplest scenario is to 
have a spin-singlet spinon pair condensate which reduces 
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U(l) to Z 2 at low energy by the Higgs mechanism—. Dif¬ 
ferent Z 2 spin liquid phases within this formalism may 
be described by inequivalent mean-field Hamiltonians, 
which can be classified by the projective symmetry group 
methocP^H. The mean-field states can be converted to 
physical spin states by Gutzwiller projection and can in 
principle be studied numerically by the variational Monte 
Carlo method^. 

This paper is outlined as follows. In Section |T[ we re¬ 
port the results of PSG analysis for square lattice and 
the mean-field representation of Z 2 spin liquids which are 
potentially relevant to the J 1 -J 2 model. We then discuss 
the properties of these Z 2 spin liquids under mean-field 
approximation in Section [TT| Section [IV| contains further 
discussion and a summary of results. Technical and nu¬ 
merical details are put in the appendices. 

I. PROJECTIVE SYMMETRY GROUP 
ANALYSIS OF SQUARE LATTICE SCHWINGER 
BOSON STATES 

The projective symmetry group analysis is based on 
a specific mean-field theory. A mean-field treatment of 
square lattice Heisenberg model and its PSG analysis us¬ 
ing fermionic spinons have been studied by Wen^S. In 
this paper we will study the spin liquid states on square 
lattice using the Schwinger boson mean-field theory. 

With the help of the Schwinger boson representation 
of the spin operator © and SU(2) completeness relation 
<j aa ' ■ cpp' = ^Sap'da'p — 5 aa i 5 p/ 3 >, the Heisenberg inter¬ 
action can be rewritten as 

Si • S, = B^Bij - 44 (4) 

where the pairing term A^ = (1/2) e aa 'bi a bja' and 

OL,Ot' 

hopping term B^ = (1/2) b\ a bj a are both SU(2) in- 

a 

variant bond operators. 

After a Hubbard-Stratonovich transformation, the 
quartic terms are decoupled and a mean-field Hamilto¬ 
nian is obtained: 

Hmf = 4 ' : A '.i A 'J + B iAo + Hc 0 

ij ( 5 ) 

+ 'y ^ (i-Ajji — \Bij\ )/Jxj — al y — ^0) 

ij i 

where complex numbers = — Aji and B,j = B* t are 
called the mean-field ansatz and the real Lagrangian mul¬ 
tiplier Hi is introduced to enforce the average boson num¬ 
ber (fii) = k = 2 S condition on every site. This mean- 
field treatment can become a controlled approximation in 
the large-V limiPSHID we will however not pursue this 
direction. 

Minimizing the variational energy with respect to the 
ansatz A, i; , B t j and fii yields the self-consistent equa¬ 
tions: 

) (Aij) = Aij/Jij , ( Bij) = Bij/Jij. (6) 


In the mean-field level, we assume that the Lagrangian 
multipulier Hi is independent of site i and the bonds A l:j 
and Bij related by symmetry operations have the same 
amplitude. 

This mean-field theory has an emergent U(l) gauge 
symmetry, namely that a local gauge transformation 

b ja -a e**Wb ja , a =t, 4-> (7a) 

A i:j —A e^+^Aij, (7b) 

—A e <[-*«+*(j)]B y , (7 C ) 

will leave the physical observables unaffected. 

The emergent gauge symmetry makes the symmetry of 
the mean-field ansatz not manifest: a mean-field ansatz 
after symmetry operation might seem different from the 
original one, but if they are connected by a gauge trans¬ 
formation the two ansatz still describe the same physical 
state and hence should be considered identical. 

In order to solve this problem, Wen and collaborators 
suggest that one should use the projective representation 
of the space group to classify different kinds of mean-field 
ansatz. 

In the projective representation, every symmetry oper¬ 
ation X is accompanied by a U(l) gauge transformation 

Gx, 

G x X(b ja ) -A < ,0 ' xu b SU)n . (8) 

The mean-field ansatz should be invariant under the 
combined operation G x X instead of symmetry opera¬ 
tion X alone. The collection of all the combined oper¬ 
ations which leave the mean-field ansatz invariant form 
the projective symmetry group (PSG). 

Different PSGs characterize different kinds of spin liq¬ 
uid states all sharing the same symmetry. Under certain 
gauge, the PSG can be fully determined by the commu¬ 
tative relations of the generators of the symmetry group. 

The PSG defined through the algebraic relations of the 
group generators is called the algebraic PSG since there 
might be PSGs that cannot be realized by any mean-field 
ansatz. 


A. PSG classification of the symmetric spin liquid 
on square lattice 

We set up a rectangular coordinate system and rep¬ 
resent the lattice site with two unit vectors e x and e y : 
r = xe x + ye y , where x, y are integers. 

The space group of the square lattice is generated by 
translation Tj along e x , translation T 2 along e y , a reflec¬ 
tion a around X-axis and the 90° rotation C 4 around the 
origin (x,y) = (0,0). 

The action of four generators on a lattice site therefore 


(fii) = K 
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FIG. 1: The coordinate system and space group generators 
T u T 2 ,C4,,a of the square lattice. 


reads (see also FIG. [I]): 


T\ 

(x,y) -t (x + i,y), 

(9a) 

t 2 

(x,y) ->• {x,y + 1 ), 

(9b) 

a 

(x,y) -t (x,-y), 

(9c) 

C4 

(x,y) -t {-y,x). 

(9d) 


These generators have the following commutative rela¬ 
tions, which completely define the space group, 


t x t 2 = t 2 t u 

( 10 a) 

Tier = crTi, 

( 10 b) 

T 2 cr = crT ' 2 -1 , 

( 10 c) 

a 2 = l, 

(lOd) 

b 

II 

b 

(lOe) 

Q 4 = 1 , 

(lOf) 

T 1 C 4 = C a T-\ 

( 10 g) 


Besides the space group generators, we can also add 
the time-reversal operator T which is commutative with 
all the space group generators. 

We have solved the algebraic PSG in the Appendix [A] 
Here we will only list the results, 


<f>T 1 

= 0, 

(lla) 

^T 2 

= PlTTX, 

(lib) 

4*17 

P47T 

= P2TTX+P37Ty + — , 

(11c) 

fat 

/ n PTK 

= Pinxy + (P2 +P 3 )^y+ ~2~, 

(lid) 

<t>T 

= ps{x + y - i)7T. 

(lie) 


where the number Pi,P2,P3,P4,P7,P8 are numbers which 
can be either 0 or 1. Therefore there are at most 64 kinds 
of PSGs. 


B. Physical realization of the PSGs on square 
lattice 

Demanding non-vanishing ansatz on certain bonds will 
impose further constraints on the algebraic PSGs. In 
Appendix[B]we have analyzed the constraints imposed on 
PSGs when demanding an arbitrary non-vanishing bond. 
Here for simplicity and with the J1-J2 model in mind, we 
will only report the results for nearest-neighbor(NN) and 
next-nearest-neighbor(NNN) bonds. 

Because of the strong antiferromagnetic Ji coupling, it 
is reasonable to assume that the nearest-neighbor pairing 
ansatz Ai is nonzero. Under this condition there are two 
classes of spin liquid states distinguished by the gauge- 
invariant flux $ in the elementary plaquette, the zero- 
flux states with $ = 0 and the 7r-fiux states with $ = tt, 
where $ is defined modulo 2 n on all plaquettes a^2l 

A ij (-A* jk )A kl (-A? i ) = \A 1 \*e i *. ( 12 ) 

In the PSG language, the quantum number of flux >I> 
is pi, with zero-flux states corresponding to pi = 0 and 
7r-flux states corresponding to p\ = 1. 

The NN pairing ansatz is invariant under a staggered 
U(l) gauge transformation, it is necessary to gap out this 
low energy U(l) gauge held by either nonzero NN hop¬ 
ping Bi or nonzero NNN pairi ng A 2 l which will reduce 
the low energy gauge held to Zj 22 * 28 [ 

The existence of nonzero Ai demands that P2 + P3 = 1 
and P 4 = p 2 - Nonzero A 2 requires that p 1 = 1 and 
P 4 + P 7 = 1- Nonzero B\ requires that p 2 = 0 and P 3 = 
1, and that B\ must be pure imaginary. Nonzero B 2 
requires that pi = 0 and B 2 is real. 

The simplest zero-flux Z2 states can be constructed 
if we demand the existence of NN pairing A\ and NN 
hopping Bi , with possible existence of NNN hopping B 2 . 
NNN pairing A 2 is however forbidden in the zero-flux 
states. 

Therefore the PSG for zero-flux states becomes 


fa 

= 4>t 2 = 0 , 

(13a) 

4^7 

= tt y, 

(13b) 

4>c 4 

, P 7 

= n V + 

(13c) 

4>t 

, u 

= {x + y- 2 ) 7r - 

(13d) 


We are still left with a Z 2 number characterizing two 
different kinds of spin liquid states. We can further 
distinguish these two states with the help of another 
gauge-invariant flux defined on all plaquettes: $9 = 
Arg [ AlJ (-A* k )(-B* kl )(-B; z )}. 

Therefore we obtain two zero-flux ( p\ = 0) spin liquid 
states: 


1 . pi = 0 , Z 2 [0,0] state with fluxes $1 = 0 and $2 = 0 
through a plaquette. A\ and B 2 are real, B\ is pure 
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imaginary. The ansatz are given by 

A( x ,y)(x+l,y) = ~^(x,y)(x,y+l) = Al, 

-®(x,y)(x+l,y) = ~^(x,y)(x,y+1 ) = B\, (14) 

-®(x,y)(x+l,y+l) = — -®(x,y)(x-l,y+l) = &2- 

2. p 7 = 1, Z 2 [0,7t] state with fluxes = 0 and $ 2 = tt 
through a plaquette. A\ and B- 2 are real, B\ is pure 
imaginary. The ansatz are given by 

A(x,y)(x+l,y) = ^4(x,y)(x,y+l) = ^1> 

-®(x,y)(x+l,y) = — -®(x,y)(x,y+l) = B\, (15) 

-B(x,y)(x+l,y+l) = ~~B(x,y)(x—l,y+l) ~ B2■ 

We can also obtain 7 r-flux (pi = 1) states by demand¬ 
ing that A-[ and A 2 are non-vanishing, with possible 
nonzero By- NNN hopping _B 2 is forbidden in 7 r-flux 
states. 

Therefore the PSG for 7 r-flux states becomes 


3. P 2 = 1, ps = 0: Z 2 [tt,tt]TZ with fluxes 4>i = n 
through the plaquette and $3 = tt through the tri¬ 
angle. Ay and A 2 are real. 

(—l) V A( x ,y)(x+l,y) = ~^4(x,y)(x,y+1) = ^1) 

( l)^ J 4(x,y)(x+l,y-t-l) ( l)^^4(x,y)(x — l,y+l) -‘4-2* 

(19) 

4. p 2 = 1, ps = 1: Z 2 \tt,tt\I with fluxes $1 = tt 
through the plaquette and $3 = tt through the tri¬ 
angle. A 1 is real and A 2 is pure imaginary. 

(— l) y 4l(x,y)(x+l,y) = ^4(a;,j/)(a:,j/-|-l) = -^l) 

(— l) !/ ^4(x,y)(x+l,y-|-l) = ( — l) i/ ^l(x,y)(x-l,y+l) = A 2 . 

( 20 ) 

II. MEAN-FIELD THEORY RESULTS 



= o, 

(16a) 

<^t 2 

= 7 nr, 

(16b) 

(j)(7 

= p 2 ttx + (1 - p 2 )ny + 

(16c) 

<t>C4, 

(1-P 2 )tt 
= irxy + ny + , 

(16d) 

4>r 

, 1 , 

= P%(x + y- -)tt. 

(16e) 


There are four 7r-flux states depending on the values of 
P 2 and ps- We can further distinguish the four states with 
gauge-invariant flux <J> 3 = Arg [Aij(-A* k )A k i(-A^)} de¬ 
fined on the triangle i —> i + x —> i + 2x —> i + x + y —> i. 

Nonzero B\ can only be realized in the p 2 = 0, ps — 1 
state since the coexistence of Ay and By requires that 
P 2 = 0 and By is pure imaginary. 

The NNN pairing A 2 can be real (corresponding to 
Ps = 0) or pure imaginary (ps = 1), which is denoted by 
the letter 1Z or X respectively. 

Therefore we obtain four 7r-flux states: 

1. P 2 = 0, ps = 0: Z 2 [tt,0]TI with fluxes flq = tt 
through the plaquette and <b 3 = 0 through the tri¬ 
angle. Ay and A 2 are real. 

( — l) 2/ ^4(x,y)(x+l,y) = J 4(x,y)(x,y+1) = 

( — l) 3/ ^4.(x,y)(x+l,y-|-l) = ( — I) 1 * ^4(x,y)(x—1,-y+l) = ^2- 

(17) 


In this section we will study the properties of the two 
zero-flux states and four 7r-flux states in the mean-field 
level. 

We will treat the average boson density k and a = 
J 2 /J 1 as variational parameters and obtain the phase 
diagram of the six states with respect to them. 

Since spinons can only be created in pairs, the phys¬ 
ical spin excitation spectrum should be the two-spinon 
continuum spectrum. We compute the lower boundary 
of the two-spinon spectrum at every given total momen¬ 
tum k for the six spin liquid states, through which we 
can distinguish different kinds of spin liquid states. 

Another measurable quantity is the static spin struc¬ 
ture factor. It can be calculated using the formula in the 
mean-field level: 

<sw = ^-E< s ‘- s i) eik ' (ri_Pi) ’ ( 21 ) 

JVsite 

t -,3 

where N s - Ite is the number of sites. 

The two-spinon spin excitation spectrum and static 
spin structure factor can in principle be measured ex¬ 
perimentally by (inelastic) neutron scattering, and nu¬ 
merically by measuring spin-spin correlation functions. 


A. Zero-flux state 


2. P 2 = 0, ps = 1: Z 2 [tt, 0]I with fluxes $1 = 7 r 
through the plaquette and $3 = 0 through the tri¬ 
angle. Ay is real and A 2 is pure imaginary. In all 
four 7r-flux states By can be nonzero only in this 
one, and must be pure imaginary. 

( — l) !/ ^4(x,y)(x+l,y) = J 4(x,y)(x,y+ 1 ) = Ay, 

( — l) y ^4(x,y)(x+l,y+l) = (~ 1) W ^4(x,-y)(x—1,-y+l) = ^ 2 , 

-®(x,y)(x+l,y) = — -®(x,y)->(x,y4» 1) = By- 

(18) 


1. ^2 [ 0 , 0 ] state 


After Fourier transformation b r = -j= e tkr & k , the 

s r 

mean-field Hamiltonian becomes 
Hmf = ^ 4>] < T) k 4' k 

k 

2 |A 1 | 2 - 2 |B 1 | 2 2 |H 2 | 2 


■ N s [fi + pn ■ 


J2 


( 22 ) 
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where we have used the Nambu spinor d'k = (£>k-t, ^-k^) T 
and the 2 x 2 matrix 


Ac = (-2 B 2 /i - n)l ~ A\f-<j y - lm(Bi)f_a z , (23) 


where f\ = sin(fc a; )sin(A:y), f± = sin(fc x ) ± sin(fc y ), 1 is 
the 2x2 identity matrix, a VtZ are Pauli matrices. 

After a Bogoliubov transformation, the mean-field 
Hamiltonian can be diagonalized to yield 


Hmf = ^2 Wk(7kt7kt + 7k|7k; + 1) 

k 

2|Ai| 2 - 2|A| 2 2|J 


+ N a [n + fin + 


Jl J 2 r 

The dispersion relation has two branches u±, 


(24) 


w k ± = \J (2£? 2 /i + /i) 2 — AJ/ 2 ± ImA/_. (25) 


The minima of dispersion are located at ±Q where 
Q = (7 t/2 , — 7 t/2 ), at which points there is an energy 
splitting between the two branches which is proportional 
to |A|- 

The self-consistent equations are 


4J3 2 /J 2 = 


f I 

<9w k + 

Jbz 2 

dn dfi 

[ I 

'<9w k + <9w k -’ 

Jbz 2 

dA x dA x 

f 1 

'<9w k + 3w k -’ 

Jbz 2 

A 0|A|j 

f I 

'9w k+ <9w k _j 

Jbz 2 

d B 2 dB 2 


d 2 fc, 


d 2 fc, 


d 2 fc, 


d 2 fc. 


(26a) 

(26b) 

(26c) 

(26d) 


It can be seen from the dispersion relation Eq. (25) that 


the ansatz Bi does not enter the self-consistent equations 
in the mean-field level. A closer examination also shows 
that B i does not affect the Bogoliubov transformation 
used to diagonalize the Hamiltonian. But B\ does have 
effect on the low-energy effective field theory and mag¬ 
netic ordered states as we shall see later. 


2. Z 2 [0, 7rj state 


The mean-field Hamiltonian after Fourier transforma¬ 
tion becomes: 


Hmf — ^2 ^kAc^k 


+ Ng [fi + fXK + 


2 |A 1 | 2 - 2 |H 1 | 

Ji 


21 Bo 


(27) 


J 2 


where 


A = (-2A/i - H)1 - A ± f + a y - Im(A)/-^, (28) 

The spinon dispersion relations are Wk± = 

Jimfi + M) 2 - A\f\ ± Im(A)/-. 


K 



FIG. 2: The mean-field phase diagram for the zero-flux states 
as a function of ratio a = J 2 /J 1 and the average boson density 
k. The two kinds of zero-flux states cannot be distinguished 
by their energetics since the NN hopping ansatz B 1 does not 
enter the mean-field equations. The curves of critical k c are 
plotted. This state has lower energy than the 7r-flux states 
in the region where a is relatively low compared to unity. 
Above the critical value of k, magnetic order will develop at 
several k points where the spinon dispersion becomes zero. 
The magnetic ordered states for the two kinds of zero-flux 
states are different. The magnetic order obtained from the 
^[0,0] state after condensation of boson is the canted Neel 
order, and the magnetic order from the Z 2 [0,7t] state is the 
Neel order. 



FIG. 3: Static structure factor for the zero-flux state. The 
spin structure factors for the two zero-flux states are qualita¬ 
tively similar so here we only show the static structure factor 
for Z 21 0, 0] state at k = 0.3. Axes are k x and k y in dimension¬ 
less units and the static structure factor is calculated using 
Eq. ( |21| ). The global maxima are located at two wavevec- 
tors: (0, 0) and When the phase transition from the 

zero-flux Z 2 spin liquid states to the magnetic ordered states 
happens, the Bragg peak will be located at these two wavevec- 
tors, which represents the Neel order. 


The minima of dispersion are located at ±Q where 
Q = (7t/2,7t/2). 

The self-consistent equations are Eq. (26a)-(26d). The 
bond Bi serves as the Higgs field to break the U(l) gauge 
symmetry down to Z 2 . Since B\ does not enter the self- 
consistent equations, the two kinds of zero-flux states can 
not be distinguished by their energy at least in the mean- 
field level. 

We have obtained the phase diagram of the zero-flux 
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state in FIG. [2j treating average boson density n and 
a = J 2 /J 1 as variational parameters. The curve of crit¬ 
ical boson density n c is plotted for the zero flux state, 
above which the energy gap closes and magnetic order is 
developed. 


B. Magnetic order from zero-flux state 


In the Schwinger boson formalism, magnetic order is 
obtained via the condensation of bosons. When boson 
density k exceeds critical value k c . spinon dispersion will 
become zero at several Q points. Bosons will condense 
at these Q points and develop magnetic orders. 

We have numerically computed the static spin struc¬ 
ture factor for the two kinds of zero-flux states and find 
that there is no qualitative distinction between the two 
cases. Therefore we only show the static structure factor 
of the Z 2 [0,0] state at a relatively low k = 0.3 in Fig. [3] 
We know from Fig.[3]that the global maxima of the static 
structure factor are located at (0,0) and ( 7 r, 7 r), which in¬ 
dicates that the spin liquid state is adjacent to the Neel 
ordered state when the phase transition from the zero- 
flux Z 2 spin liquid states to the magnetic ordered states 
happens. 

The analytical analysis of magnetic order from Z 2 [ 0,0] 
and Z 2 [ 0 , 7 r] states are given in the Appendix [C| Here 
we will report only the final results. For convenience of 
later discussion, we devide the square lattice into two 
sublattices, sublattice u for (—1 ) x+v = 1 and sublattice 
v for (-1 ) x+ y = - 1 . 

The magnetic ordered state obtained from the Z 2 [ 0,0] 
state is a non-collinear order, given by 


<S,r,) . 


1 - (-1 ) x+ v 


(29) 


where n UjV are two vectors of the same length, \n u \ = 
n v |. The angle between n u and n v are in general 
less than 180° unless B 1 = 0. This represents the 
“canted Neel order”, which is the classical ground state 
for nearest-neighbor Heisenberg model under a small uni¬ 
form magnetic field. 

Aside from the spin (magnetic dipole) order parameter 
n UtV , this canted Neel order also has another vector spin 
chirality order parameter 


C i:i =v ij {S i x Sj), (30) 


where Vjj = ±1 on the nearest neighbor bond (■ ij ), and 
takes value + 1 (— 1 ) if site i belongs to sublattice u(v). 

The magnetic order obtained from the Z 2 [0,tt] state is 
the collinear Neel order, (S(r)) = (— 1 ) x+v n. But this 
state is not the conventional Neel order and will also 
acquire a nonzero ground state expectation value of the 
vector spin chirality order parameter (S, x Sj) due to 
quantum fluctuations. 

In the effective field theory, the transition from Z 2 
spin liquid state to the magnetic ordered state occurs 


through the condensation of the charge -1 [in terms of 
the emergent U(l) gauge field] spinon field z a . B\ serves 
as a charge- “-2” Higgs field. The physical observables 
are charge-“0” combinations of the spinon held z a . The 
three possible bilinear combinations are as follows: 

ni = z'Sz, n 2 = Re[f?i z T ia v az\, n 3 = Im[Hi z T ia v az], 

(31) 

which are orthogonal to each other. 

It is easy to see that Hi can be identified with Neel or¬ 
der parameter. And since the vector spin chirality order 
parameter is even under time-reversal transformation, we 
can identify it with Re[i?i z T icr y az\. The order parame¬ 
ter is in fact the triad ^ 1 , 2 , 3 , even if the expectation val¬ 
ues of spins form the collinear Neel order. This kind of 
magnetic order has been obtained before on honeycomb 
lattice in proximity to certain spin liquid state 31 . 


C. Critical field theory for transition from 
zero-flux spin liquid states to magnetic order 


When k exceeds the critical value k c , there will be a 
phase transition from the Z 2 spin liquid state into the 
magnetic ordered state. Close to the phase transition 
point a continuous field theory can be derived. 

The low-energy effective Lagrangian for the transition 
from Z 2 [ 0,0] state to the canted Neel order reads 


£= d?v{ -————— —Im{Bi)uj* a —uj a 

J (—fi + 2A\ + 2B2)cl z clt 


+ (Ai - 2 B 2 + 

+ ( 


4|Ri 


—[i + 2Ai + 2B 2 
—2 fj, — 4Ai -f- 4lB 2 8|n>i 


)<9 r w* • d r uj a 


8|Bi' 2 


(—n + 2 Ai + 2 B 2 )a 2 


)w> Q }, 

(32) 


where ui is related to the spinon field z and is defined in 
Eq. (D10). Summation over repeated indices is implied 
hereafter. 

The coefficient of the mass term u*w a will change sign 
upon approaching the critical point, indicating that the 
critical chemical potential is /x —>■ 2B 2 — ^/4 Af + 4|£>i| 2 . 
The low-energy effective field theory now flows to a fixed 
point where space and time scales differently with dy¬ 
namical critical exponent z = 2 , which is different from 
previous theories of such phase transitions?^. By power 
counting one can find that this theory approaches the 
upper critical dimension d = 4 where the physics is 
controll by a Gaussian fixed point and the fluctuation 
around the mean-field theory is negligible. The scal¬ 
ing properties of correlation functions can be determined 
by naive power counting since in the upper critical di¬ 
mension, the anomalous dimension correction approaches 
zero in the 4 — e analysis when e —> 0. For example, 
(S(0) • S(r)) oc where 77 is 4 plus correction which 

is proportional to e and hence is negligible. 
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The low-energy effective Lagrangian for the transition 
from Z 2 [ 0 , 7 r] state to the Neel ordered state is 

£ = y d 2 r{\D T z\ 2 + c 2 \D T z\ 2 + m 2 \z\ 2 }. (33) 

The couplings between Higgs and spinon fields have 
at least two spatial derivatives and one time derivative, 
which is irrelevant by naive power counting. Anomalous 
dimension will not change this result. Upon approach¬ 
ing the critical point, the mass term will disappear and 
this theory has an enlarged 0(4) symmetry. The scal¬ 
ing properties of the field theory at 0(4) critical point 
have been studied both analytically and numerically. The 
scaling behaviour of the correlation functions is therefore 
known. As an example, the spin-spin correlation func¬ 
tion has a relatively large anomalous dimension, which 
behaves as (S(0) • S(r)) oc |r| —77 , where ?y is numerically 
determined as 77 = 1.373(3) in contrast to the result for 
Z 2 [ 0,0] stat^ 33 : 

The transformation rules for spinon fields (ui or z) un¬ 
der space-group symmetry can also be readily deduced 
from the PSG and are listed in TABLE|l| It is easy to ver¬ 
ify that the form of Lagrangian is invariant under these 
transformation rules. 


Z 2 [0, 0] state 

Z 2 [ 0,7r] state 

Symmetry Operation 

uj —> —iuj 

2 — > —ia v z* 

Ti 

uj —y iuj 

z —¥ ia y z * 

t 2 

uj —y uj 

z z 

c 4 

LJ —»■ UJ 

z —y z 

a 

uj —y uj 

z — > a v z* 

T 


in which, 


Ac = 


[i ■ 1 A-, P\ + A 2 P 2 

—AiPi — A^P 2 y 1 


And we have also used the Nambu spinor 


= 


^ frukt \ 
by kt 
h t 

"u-kj. 


and two 2x2 antihermitian matrices 


1 / 2isin(fc x ) — l + e 2zfcl A 

2 yl — e~ 2tk v —2isin(k x )J ’ 


(36) 


P'2 


1 ( 0 

2 0 )' 


(37) 


- g^^a: g ikx | g ^(^a: I ~‘ 2k y') gd^a; <2k y'} 

The Hamiltonian can be diagonized by a Bogoliubov 
transformation using a SU(2,2) matrix to yield : 


Hmf = J2 w k(7Lt7akt + + 1) 

k;a=M,u (38) 

+ At/ 7 + /tre + 2|Ai|“/ J\ + 2\A 2 \ 2 / J 2 \. 

In this state A 2 is real, the dispersion relation is there¬ 
fore fourfold degenerated: 


TABLE I: Transformation rules of the spinon fields under 
PSG. 


D. 7r-flux state 


w k = \jlP ~ A\fi - 4 A 2 / 2 , (39) 

where /1 = sin 2 (fc x ) +sin 2 (fc y ), / 2 = sin 2 (/c x )sin 2 (fcj / ). 

The minima of spinon dispersion are located at 
(■ k x ,k y ) = ±(7r/2, 7 t/2). 

The mean-field self-consistent equations are as follows: 


1. Z 2 [n, 0 ] 1 Z state 

There are two sites in a unit cell of 7r-flux ansatz dis¬ 
tinguished by (— l) y = ±1 which can be labeled by u and 
v respectively. The unit cells are labeled by integers x 
and y , where the u(v) site in the unit cell at ( x , y) are at 
position (x , y) = ( x , 2 y) [{x, y) = (x, 2 y + 1)]. 

After Fourier transformation 

b, . = 1 V^p-hfcxX+2 k v y) b 734') 

u,v)rot fAf~ / j u {u,v)K.or> 

V S x,y 

where k X:V = k ■ e XiV , the mean-held Hamiltonian be¬ 
comes: 

Hmf = 'y ' fE , j c ^k^ r k±-N s [y+ii k+2\A\ | 2 / Ji±2| A 2 1 2 / J 2 ], 

( 35 ) 


1 + re = — [ d 2 fc, (40a) 

Jbz d n 

4A 1 /J 1 = ~ [ f^d 2 fc, (40b) 

Jbz oAx 

AA 2 /J 2 = - [ d 2 fc. (40c) 

JBZ U A 2 

Solving the self-consistent equations we have obtained 
a mean-held phase diagram for this state as shown in 
FIG. [t] This state has re c ss 0.81 but exist in the regime 
a = J 2 /Ji ~ 0.75 ~ 1 far away from the physically inter¬ 
esting regime a ~ 0.5. Therefore this state is unlikely to 
be the physical ground state for the square lattice J\ — J 2 
Heisenberg model found in the numerical methods. The 
lower edge of the two-spinon spectrum is displayed in 
Fig. [5] for the Z 2 [tt,0]1Z state at re = 0.5, which could be 
used as a probe for this state in numerical simulations. 














K 


2. Z 2 [tt,Q]Z state 



FIG. 4: The mean-field phase diagram for the Z2[tt, 0]TZ n- 
flux state as a function of ratio a = J 2 /J 1 and the average 
boson density k. This state has a critical value of ft fss 0.81, 
and occupy a finite area in the regime a = J 2 /J 1 ~ 0.75 ~ 1. 
Upon decreasing the value of a, the Z 2 spin liquid state will 
become a U(l) spin liquid state with only non-zero nearest- 
bond A\ through a first-order phase transition. The U(l) spin 
liquid state with only A\ has a critical value of k ss 0.81. Upon 
increasing the value of a, the Z 2 spin liquid state will become 
a U(l) spin liquid state with only non-zero A 2 through a first- 
order phase transition. The U(l) spin liquid state with only 
A 2 can be regarded as two copies of decoupled 7r-flux states 
on the square lattice and has a critical value of k ss 0.40. 



This state has nonzero B\, so after Fourier transfor¬ 
mation, the Hamiltonian becomes 


Hmf - 

k 


/ (i-l + BiP 3 
y-A^ - A* 2 P 2 


At p i + A 2 Po\ 
ft ■ 1 - B ± P 3 j 

(41) 


FIG. 5: The lower edge of the two spinon spectrum for the j n g^ate A 2 is pure imaginary, we have two twofold- 

Z 2 [tv,0]TZ state at k = 0.81. degenerate dispersion : 


^k± = \jft 2 - (Aj | B^)h - 4|Z| 2 h ± 2^/(| Btf + 4(A 2 - |5 1 | 2 )|A 2 p/ 2 )/i, 


(42) 


where /, = sin 2 (/c a; )-|-sin 2 (fcj / ) ) f 2 = sin 2 (fc x )sin 2 (fcj,) and 


1 / 2isin(fc x ) — 1 + e 2iky \ 

2 yi — e - 2lk y —2isin(fc x ) J 


(43) 


Minima of the spinon dispersion are at ( k x ,k y ) = 
±(7r/2,7r/2). 


The self-consistent equations are as follows: 


1 + k = — 



AA\j J\ — — I 
Jbz 


415,1/J, 

4|H 2 |/J 2 




(44a) 


(44b) 

V W k+ , 

(44c) 

l[ 9lJ k + , d ^-]P h 
2 [ 9|h 2 | + d\A 2 y d k - 

(44d) 
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The self-consistent equations are Eq. (40a)-(40c). The 
mean-held phase diagram for the Z 2 [ir, OJI state is ob¬ 
tained as shown in FIG. [6j This state has a relatively 
low k c . The lower edge of the two-spinon spectrum is 
displayed in Fig. [7] for the Z 2 [n, 0]Z state at n = 0.6. 


where = c^ kx T e ^ kx -t- e 

The dispersion relation is fourfold degenerate: 


w = \Jn 2 - A\f Y - AAjfo , (48) 


K 



FIG. 6: The mean-field phase diagram for the Z 2 \f : Q]T 7r-flux 
state as a function of ratio a = J 2 /J 1 and the average boson 
density k. This state has a critical value of k « 0.60, and 
occupy a finite area in the regime a>0.10. 



FIG. 7: The lower edge of the two-spinon spectrum for the 
^ 2 [ 70 , 0]I state at k = 0.60. 


where /1 = sin 2 (fc 2; ) + sin 2 (£: y ), / 2 = cos 2 (fc a; )cos 2 (fcj / ). 

We denote |A 2 |/Ai = tan(0). When f?<arctan(l/-\/2), 
the minima are located at ±(7 t/2 ,7 t/2 ), and when 
0 >arctan(l/-\/2), the minima jump to ( k x ,k y ) = (0,0) 
and (7T, 0). 

The curve of critical k is obtained via the self- 
consistent equation as shown in FIG. [8] This state has a 
relatively high value of n c « 1.36 even above the physical 
value k = 1. Close to the physically interesting regime 
where k = 1 and a = J 2 /J 1 ~ 0.5, this state occupy a 
finite area of phase space, and hence might be a promis¬ 
ing candidate for the numerically found Z 2 spin liquid 
state. But a closer study shows that this state is ener¬ 
getically more unfavorable than the Z 2 [n, 7r]I state in the 
physically interesting regime. The lower edge of the two- 
spinon spectrum is displayed in Fig. [9] for the Z 2 [tt,tt]71 
state at k = 1.36. 


K 



FIG. 8: The mean-field phase diagram for the Z 2 [n,n]lZ n- 
flux state as a function of ratio a = J 2 /J 1 and the average 
boson density k. This state has a critical value of k « 1.36, 
and occupy a finite area in the regime a = J 2 /J 1 ~ 0.5 ~ 1.9. 


3. Z 2 [n, 7r ]TZ state 


The Hamiltionian after Fourier transformation is now 
(up to a constant): 


H mf = 

k 


( Ai-1 

V-^rPi - a* 2 p 2 


A 1 P 1 + A 2 P 2 
!-<■ ■ 1 , 


And we have used two 2x2 matrices: 


*k- 

(45) 


1 / 2ism(k x ) 1 — e 2lky \ 

2 1 + e~ 2lky —2ism(k x )J 


(46) 



Pi = 


1 

2 




(47) 


FIG. 9: The lower edge of the two-spinon spectrum for the 
Z 2 \k, tv]1Z state at k = 1.36. 
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4- Z 2 [ir, n]T state 


The mean-field Hamiltonian is Eq. (45). 

In this state, A 2 is imaginary, so the dispersion rela¬ 
tions are twofold degenerated: 


W± = yjv 2 - A\h - 4|H 2 |2/ 2 ± (49) 


where /i = sin 2 (fc x ) + sin 2 (fe y ), / 2 = cos 2 (fc a; )cos 2 (fcj ; ). 

Define |A 2 |/A! = tan(0). When 0<arctan(l/(2v / 2)), 
the minima are at ( k x ,k y ) = ±{n/2,n/2). When 
0>arctan(^j), the minima will move to incommensu¬ 
rate wave vectors (k x ,k v ) = ±(q,q), where sin (q) = 
l/(2 v / 2tan(6»)). 


The self-consistent equations are Eq. (44a)-(44d I and 
the mean-field phase diagram is obtained as shown in 
FIG. HO] 

This state has a high n c k, 1.62 even above the phys¬ 
ical value k — 2S = 1 and occupy a finite area of phase 
diagram close to the physically interesting regime n ss 1 
and a « 0.5. Although this state is energetically more 
unfavorable than the zero-flux states in this regime, we 
propose that by adding ring-exchange term in the Hamil¬ 
tonian may favor the 7r-flux spin liquid state. The lower 
edge of the two-spinon spectrum is displayed in Fig. ED 
for the Z 2 [7r,7r]I state at k = 1.62. 



FIG. 11: The lower edge of the two-spinon spectrum for the 
Z 2 [n,n]X state at k = 1.62. 



K 



FIG. 10: The mean-field phase diagram for the Z2[7r,7r]iT n- 
flux state as a function of ratio a = J 2 /J 1 and the average 
boson density k. This state has a critical value of k s=s 1.62, 
and occupy a finite area in the regime 0.55 < a < 1.27. 


E. Magnetic order from 7r-flux states 


We have computed the static spin structure factor for 
the four kinds of 7r-flux states in the mean-field level and 
find that there is no qualitative distinction between these 
four cases. Therefore we only show the static structure 
factor of the Z 2 [7r, 0]7£ state at a relatively low n = 0.2 
in FIG. 12 From Fig. [12] that the global maxima of the 
static structure factor are located at ±(7r, 0) and ±(0, 7r), 
which indicates that the magnetic ordered state adjacent 
to the 7r-fhix state is different from the Neel order. 


FIG. 12: Spin structure factor for the 7r-flux state. The spin 
structure factors for the four 7r-flux states are qualitatively 
the same so here we only show the static structure factor for 
^ 2 [tt, 0]77 state at k = 0.2. Axes are in dimensionless units 
k x and k y and the static structure factor is calculated using 
Eq. |21| The global maxima are located at four wave vectors 
±(-7r,0) and ±(0,7r). It is expected that when the continuum 
phase transition from the 7r-flux Z 2 spin liquid states to the 
magnetic ordered states happens, the Bragg peak will located 
at these four wave vectors, which represents a new kind of 
magnetic order different from the Neel order. 


We have also analytically obtained the magnetic or¬ 
dered state starting from the mean-field Hamiltonian us¬ 
ing the method of boson condensation. For simplicity, we 
shall only consider a 7r-flux state which only has NN bond 
A-[. The mean-field Hamiltonian after Fourier transfor¬ 
mation is 


HMF = ^ k 

k 


( M-l. 

y - A 1 Pi, 


AiPi\ 

1 J 


*k, 


where 


1 / 2*sin(fca;), -l + e^A 

2 \ \ — e~ 2lkv , —2*sin(fc x )y 


(50) 


(51) 


The critical spinon vector occurs at two inequivalent k 
points ±Q, where Q = (tt/2, 7t/2). 
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Under the constraint that the density of condensate on 
every lattice site is uniform, we can work out the pattern 
of spinon condensation. 

We then find a set of four-sublattice ordered states 
consistent with a subset of the classical ground state for 
J 2 /J 1 = 1/2 Heisenberg model. 

The magnetic order is (S(r)) = m 1 (—I)" 1 + m 2 (— l) v + 
m 3 (—l) x+y , where 1 x 11 , 2,3 are three orthogonal vectors 
and 

m\+m\+m\ = m 2 . (52) 

Technical details of the calculation can be found in the 
Appendix |C 3| 

III. DUALITY BETWEEN SCHWINGER 
BOSON SPIN LIQUID STATES AND 
ABRIKOSOV FERMION SPIN LIQUID STATES 

In the parton constructions of the Heisenberg model, 
the physical spin operator can be decomposed into 
spin-1/2 partons that can be either bosonic (Schwinger 
bosons) or fermionic (Abrikosov fermions). Whether the 
two seemingly distinct approaches are equivalent or not 
remains a long-standing puzzle. In this section we closely 
follow the work of Hermele et and Lu et and 
deduce the correspondence between the PSGs of the 
Schwinger boson representation and Abrikosov fermion 
representation. 

We know clearly from the PSG classification of spin 
liquid states that topological order alone is not enough 
to fully characterize the different phases of the spin liquid 
states when symmetry is also presented 331 . Actually, the 
interplay of topological order and symmetry will yield a 
richer stru ctur e called the “symmetry enriched topolog¬ 
ical order’I 2 ®!! 33 3 ®2l0 i n different symmetry enriched 
topological phases, anyon excitations will not only have 
fractional charges and fractional statistics, but also carry 
fractional symmetry quantum numbers. To be more con¬ 
crete, the symmetry operation acts on anyons projec- 
tively and when an anyon returns to its original posi¬ 
tion after a series of symmetry operation, it may gain 
a nontrivial phase factor due to the gauge structure of 
the theory. All the fractionalization of symmetry quan¬ 
tum numbers of a certain type of anyon put together 
defines a fractionalization class for this type of anyon. In 
the case of Schwinger boson spin liquid states, the frac¬ 
tional symmetry quantum numbers pi (i = 1, 2,3,4, 7, 8) 
fully determine the fractionalization class for the bosonic 
spinon. Once the fractionalization classes for each type 
of anyon excitations of a topological ordered state is de¬ 
termined, we can specify a symmetry class for this kind 
of symmetry enriched topological phase. 

The fusion rules between different types of anyons have 
certain constraints on the fractionalization classes in the 
same symmetry class. Therefore we could use the com¬ 
patibility condition to determine the correspondence be¬ 
tween the fermionic PSGs and bosonic PSGs. 


In the present case, three kinds of topological excita¬ 
tions, bosonic spinon b , fermionic spinon / and the vison 
v, obey the following fusion rules. 

bx f = v, bxv = /, f xv = b, 

, , r f -1 (^ 3 ) 

bxb=jxf=vxv= 1. 

The fusion rules have a strong constraints on the sym¬ 
metry fractionalization between the three kinds of topo¬ 
logical excitations. Therefore in principle we could obtain 
the vison PSG from the knowledge of the fermionic PSG 
and the bosonic PSG due to the fusion rule b X v = f. 
The additional phase e l ^ f picked up by a fermion can be 
obtained from the phase of a boson, the phase of 
a vison, and possibly a twist factor e 1 ^ from the mutual 
statistics of bosons and visons, therefore we have 

= e i4>v . e i<tn, . e i4>t _ ( 54 ) 


Algebraic Identities bosonic b a 

fermionic f a 

vison v = b x f 

T 2 - 1 TiT 2 T 1 - 1 

(-ir 

Pxy 

-1 


(_f)P2 

flxpy 

-1 

a~ 1 T2aT2 

(_1)P3 

T]xpx 

1 

a 2 

(_1)P4 

v* 

1 

a^CiaCi 

(-ir 

VcrC4 

1 

ci 

1 

VCi 

-1 

C^'T\C,T 2 

1 

-1 

-1 

C^ 1 T 2 C i T~ 1 

(_1)P2+P3 

xpx TJxpy 

1 

C^T^CaT 

(_1)P7 

VCiT 

1 

3 

1 

HI 

1 

S 

(_1)P4 

VctT 

1 

T 1 - 1 T“ 1 TiT 

(_1)PS 

Vt 

1 

T 2 - 1 T“ 1 T 2 T 

(_1)PS 

Vt 

1 

T 2 

-1 

-1 

1 


TABLE II: Correspondence between bosonic and fermionic 
PSGs. The bosonic PSGs are labeled by six integers pt = 
0,1 (i = 1, 2, 3,4, 7, 8 ), while fermionic PSGs are labeled 
by nine integers {Vxy,Vxpx,ilxpy,Va,'nc 4 ,'n^c i ,Vc 4 T,ri>TT,'nt) 
where 77 = ±1. The vison PSG are fully determined as shown 
in AppendixjE] The identity C^ 1 T\CaT 2 is not gauge invari¬ 
ant as we can always tune the relative phase between the 
gauge transformation G Tj and Gt 2 , therefore its PSGs can 
always be fixed as shown above. The vison PSGs can be ob¬ 
tained from the product of fermionic PSGs and bosonic PSGs 
except the two cases Cf and a 2 where an extra factor of (— 1 ) 
should be taken into consideration as discussed in the main 
text. 

The Abrikosov fermion construction of the Heisenberg 
model and its PSG study is summarized in detail in 
Wen’s paper 27 . Although the solutions of the fermionic 
PSG are under certain gauge, the universal data is en¬ 
coded in the fractionalization classes as in the case of 
Schwinger boson, and therefore we obtain the fractional¬ 
ization classes of the Abrikosov fermion as listed in TA¬ 
BLE HD 

In the mean-field level, visons can be considered as 
point-like excitation located on the center of a plaque- 
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tte. Due to the Z2 nature, two visons will annihi¬ 
late each other, therefore we shall use v to denote vi- 
son creation or annihilation operator. The dynamics 
of the vison can thus be described by a Ising gauge 
theorj@23HH anc j i ts sy mmetry fractionalization can be to¬ 
tally determine d 43 * 44 ^. The vison PSGs are determined as 
listed in TABLE mi The detailed calculation is summa¬ 
rized in Appendix [E] 

In most cases, the twist factor from the mutual statis¬ 
tics of bosons and visons is trivial, but there are three ex¬ 
ceptions where the twist factor are nontrivial (e ^ = — 1 ), 
discussed in Ref. 13.11341391 In the following we sketch 
their arguments for these nontrivial twist phase factors. 

The first nontrivi al cas e is C 4 , namely fourfold rotation 
acted by four times 33 ®!. Consider the bound state of a 
vison and a bosonic spinon, which fuse into a fermionic 
spinon. When the bound state is rotated around a loop 
under the operation of C 4 , the vison will encycle the bo¬ 
son once and therefore pick up a Berry phase of 7 T 33 I 34 I 43 [ 

The second case is cr 2 , namely mirror reflection acted 
by two times, as discussed by Lu et aZ.®land Qi et al. 39 . 
Consider a system consisting of two fermions /1 and / 2 
that can be related by a mirror reflection cr. Then we 
have ct / ict -1 = e 1 ^ 1 f2, cr/ 2 cr _1 = e 1 ^ 2 f±. Therefore act¬ 
ing a twice on fi will yield a factor of = e 1 ^ 1 • e 1 ^ 2 . 
Alternatively, acting <7 once on the pair /1 and / 2 will 
yield a factor of —e 1 ^ 1 ■ e^ 2 , in which the minus sign 
originates from the exchange the two fermions under a. 

The same argument works for the boson pair and vison 
pair, but since there is no additional statistical sign as in 
the case of fermion, we notice that the extra phase factor 
after acting a twice on the bosonic spinon or vison is 
the same as the phase factor after acting er once on a 
pair of bosonic spinons or visons. If we denote the phase 
acquired by acting a twice on b(v) as (e*^”), then the 
above statement is equivalent to —e 1 ^ 1 ■ e 1 ^ 2 = e l( ^ b ■ e l< ^ v 
(since we can simply treat / as the bound state of b and 

v)■ 

Therefore we have the nontrivial fusion rule for the 
PSG of cr 2 : e^f = -e^ ■ e i<t>b . 

The third nontrivial fusion rule comes from a~ 1 TuT- 

Let’s consider a system composed of a pair of fermionic 
spinons /1 and f 2 located on the X-axis, so that they are 
symmetric under reflection cr. And we adopt another 
assumption that the two fermions can be connected by 
translation T x , which ensures the two fermions share the 
same symmetry quantum numbers. 

In this case, the system is an eigenstate of the operator 
cr, we therefore have cr/jcr -1 = e II?il /i, cr/ 2 cr _1 = e 1 ^ 1 / 2 . 
And acting cr twice on a fermion will yield a factor of 

g * 0 / = g 2 i 0 i 

Consider a pair of bosonic spinons and visons located 
on the X-axis, and we denote the eigenvalue of b(v) under 
a as e 1 ^ 2 (e 1 ^ 3 ). Acting a twice on a boson(vison) will 
yield a factor of e l<i>b = e 21 ^ 2 (e 1 ^ = e 2i< ^ 3 ). 

If we treat the fermion as the bound state of a bo¬ 
son and a vison, then from the above discussion we have 
e 2 «£i _ _g 2 i 0 2 . e 2i03_ By splitting —1 equally into the 


two fermion sectors, we have e 1 ^ 1 = ( ±i)e 11 ^ 2 ■ e 1 ^ 3 , which 
means each fermionic spinon gets an extra ±z phase in 
additional to the phase acquired by the boson and vison 
under reflection cr. But since the time-reversal transfor¬ 
mation will take a number to its complex conjugate, the 
extra phase will be (±z)*(±z) = 1 under the symmetry 
operation (<tT) 2 , which indicates that (&T) 2 has a trivial 
fusion rule. 

Consider the following algebraic identity 

(7V ) 2 = (a-'T-'oT) ■ (T ) 2 • (n) 2 , (55) 

since the fusion rule for (Ter) 2 and T 2 are both trivial, 
while the fusion rule for a 2 is nontrivial, we know that 
the fusion rule for cr~ 1 T~ 1 crT is nontrivial. 

Notice that the fourfold rotational symmetry C 4 is 
a combination of two reflection symmetry C 4 = 7A y cr, 
where P xy is defined as the reflection along y = x. There¬ 
fore from the algebraic identity 

Cz'T-'CaT = (PTyT^PxyT) ■ (a-'T-'oT), (56) 

we know that the fusion rule for C^ 1 T~ 1 CTT is trivial. 

The correspondence of fermionic PSGs and bosonic 
PSGs are listed in TABLE [Tl] From the fusion rule, 
we find Z 2 spin liquids that can be both described by 
fermionic and bosonic partons if the following conditions 
are satisfied: 

Vxy = (-1) P1 + 1 , Vxpx = (-1) P3 , Vxpy = (-l) P2 + \ 

lit = (—l) ps , f?o-r = Vo- = (-1) P4 + 1 , VaCi = VC 4 r = (-1) P7 , 
7c 4 = 1, 

(57) 


where the 77 s are ±1 which label different fractionalization 
classes of fermionic spinons (see Appendix [F] for their 
definitions). 

All the six Schwinger boson Z 2 spin liquid states 
we have discussed have fermionic counterparts. The 
Schwinger boson Z 2 [0, 0] state has two fermionic counter¬ 
parts, which are Z 2 Bt 3 t 3 t 1 t 3 state and Z2Bt 2 ^t 2 _t 1 t^_ 
state. The Schwinger boson Z 2 [ 0,7r] state, Z2[tt,0]1Z 
state and Z'T\jr, 0]Z state correspond to Z 2 Bt 3 t 3 t°t 3 
state, Z 2 At 3 t 3 t°t (1 state and Z 2 At}_t^t°t 3 state re¬ 
spectively. The Schwinger boson Z 2 [7t 7 n]lZ state cor¬ 
responds to Z 2 Ar)j(r)j(r 0 r^ and state. 

The Schwinger boson Z 2 [7r, tt]L state corresponds to 
Z 2 At 3 t 3 t°t 3 and Z 2 At 3 t)!t 0 t£ state. See TABLE III 


Their possible realization on square lattice is discussed iii 
AppendixjGj Among these nine states, the Z 2 Ar)!.T)(.T 0 r 3 
state cannot be realized by mean-field ansatz, and the 
Z 2 AtIt%t 9 t 9 _ state is a gapless spin liquid state whose 
gaplessness is protected by its PSG^L, and the remaining 
PSGs can all be realized by gapped spin liquid states. 


IV. DISCUSSIONS AND SUMMARY 

Here we discuss the relation between our results and re¬ 
lated theoretical and experimental works. The PSG clas¬ 
sification of Schwinger boson states has been used in Tao 
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Schwinger boson 

Abrikosov fermion 

(P1,P2,P3,P4,P7,P8) 

Label 

(v xy 5 7} xpx , Tjxpy , 7]t ) 


gp v 


9 T 

Label 

Perturbatively gapped? 

(0,0,1,0,0,1) 

£2 [0,0] 

(-1,-1,-1,-1) 

IT 1 

IT 1 

ir 1 

IT 3 


Yes 

IT 2 

IT 2 

IT 1 

IT 3 

ZaBririrVi 

Yes 

(0,0,1,0,1,1) 

Z2 [0, 7r] 

(-1,-1,-1,-1) 

IT 1 

IT 1 

T° 

IT 3 

Z 2 BtLtLt 0 T- 

Yes 

(1,0,1,0,1,0) 

Z 2 [ 7T, 0]7 Z 

(i,-i,-i.i) 

IT 1 

IT 1 

T° 

IT 3 

Z 2 At1t1t°tI 

Yes 

(1,0,1,0,1,1) 

Z 2 [ 7T, 0]2 

(1,-1,-1,-1) 

IT 1 

IT 1 

T° 

IT 3 

Z2AtLtLt°t1 

Yes 

(1,1,0,1,0,0) 

Z2 [tt, n]lZ 

(l,i,l,i) 

r° 

T° 

T° 

T° 

Z 2 AtZtIt°tZ 


T° 

T° 

r° 

IT 3 

Z 2 AtIt° + t 0 tI 

Yes 

(1,1,0,1,0,1) 

Z 2 [tv, 7t]X 

(i,i,1,-1) 

T° 

T° 

T° 

T° 

Z 2 AtIt 0 + t 0 t° 

No 

T° 

T° 

T° 

IT 3 

Z 2 At+t+t 0 t1 

Yes 


TABLE III: Correspondence between PSG solutions of bosonic and fermionic Z2 spin liquids. The Schwinger boson Z%[ 0 , 0 ] state 
has two fermionic counterparts, which are Z2 1 Bt\t\t 1 t^_ state and Z^BtIt^Lt 1 state. The Schwinger boson ^2[ 0 ,7r] state, 
Z 2 [iv, 0]1Z state and ^2[7r, 0 ]I state correspond to Z 2 BtLtLt°t £ state, Z2Aririr°r4. state and Z 2 At 1 t 1 t°t 1 state respectively. 
The Schwinger boson ^2 [vr , 7r] 1Z state corresponds to and Z2Ar°r+r 0 r|. state, of which the ^2Ar^r°r°r+ state 

does not have a mean-field realization 2 '. The Schwinger boson Z 2 [7T, 7 t]X state corresponds to Z 2 At+t 9 t°t 9 and Z2 At+t+t°t£ 
state. Among the eight PSGs that have mean-field realization, only the Z 2 At 9 t 9 t°t 9 is a gapless state whose gaplessness is 
protected by its PSG 51 . The remaining seven PSGs can all be realized by gapped spin liquid states. 


Li et aP^ and Yi-Zhuang You et aPl In particular Tao 
Li and coworkers studied the energetics of several kinds 
of projected Schwinger boson wave functions on square 
lattice for J\-J -2 Heisenberg modeP^. However they did 
not provide a complete projective symmetry group anal¬ 
ysis and presented numerical results for zero-flux states 
only. Yi-Zhuang You and coworkers tried to explain the 
behavior of local moments in iron-based superconductors 
by Schwinger boson spin liquid states 7 . They achieved a 
complete solution of algebraic PSG in appendix, but did 
not study all the mean-field states. As far as we know our 
paper is the first complete account on the PSG classifica¬ 
tion and mean-field realizations of Schwinger boson spin 
liquid states on square lattice, as well as the magnetic 
orders connected to Schwinger boson Z 2 spin liquids. 

We do not attempt to relate our results directly to any 
experimental or numerical evidences of spin liquids. It 
may be tempting to try to relate these Z 2 spin liquids 
with underdoped cuprate superconductors. This possi¬ 
bility has been thoroughly discussed in the review by Lee 
and coworkersP^, albeit using a different slave-particle for¬ 
malism with fermionic spinons. If the relation between 
Schwinger boson and Abrikosov fermion spin liquids can 
be firmly established, our results may also be relevant 
to cuprates. More recently an inelastic neutron scatter¬ 
ing experiment observed evidence of spinons coexisting 
with Neel order®. This type of states is beyond our cur¬ 
rent study of symmetric spin liquids without any sym¬ 
metry breaking. Classification and numerical studies of 
this type of “AFM*” states will be an interesting future 
direction. 

Finally we would like to say a few more words about 
the relation between Schwinger boson and Abrikosov 
fermion spin liquid states. 

Several groups^ 9 HHII1 have established the exact map¬ 
ping between these two formulations on extremely short- 
range(nearest-neighbor) resonating valence bond(RVB) 


states, which are limiting cases of Gutzwiller pro¬ 
jected wave functions of mean-field Schwinger boson or 
Abrikosov ferrnon states. This exact mapping however 
is not very useful on square lattice. There are two 
nearest-neighbor RVB states on square lattice preserving 
all lattice symmetry. They correspond to the zero-flux 
Schwinger boson(7r-flux Abrikosov fermion) states and 
the 7r-flux Schwinger boson(zero-flux Abrikosov fermion) 
states respectively. However the nearest-neighbor RVB 
states on square lattice are not Z 2 spin liquidspS, and 
they cannot distinguish different mean-field ansatz with 
the same flux in elementary plaquette. To represent 
Z 2 spin liquids and distinguish all PSG classes, valence 
bonds beyond nearest-neighbors will be required. We 
have checked on small lattices that the exact mapping 
for the nearest-neighbor RVB stated®®® cannot be es¬ 
tablished with long range valence bonds. 

In summary, we have studied Schwinger boson spin 
liquid states on square lattice by the projective sym¬ 
metry group analysis and Schwinger boson mean-field 
theory. Six symmetric Z 2 spin liquid states, two zero- 
flux states and four 7r-flux states, have been identified 
to be possibly relevant to the J 1 -J 2 model. The zero- 
flux states can go through continuous phase transitions 
into canted or collinear Neel order. However these mag¬ 
netic orders always show nonzero vector spin chiral¬ 
ity (Sj x Sj), even if the spin expectation values are 
collinear. The 7r-flux states are in proximity to certain 4 - 
sublattice magnetic orders including the Neel order, and 
can be favored energetically over zero-flux states with 
ring-exchange interactions^®. At mean-field level, two 
of the 7r-flux states, the Z 2 [k,'k\R. and Z 2 [7 r,7r]I states, 
may be stable against magnetic order for spin-1/2 and 
around J 2 /J 1 ~ 1/2. We have computed the bottom of 
two-spinon excitation continuum and static spin struc¬ 
ture factors for these spin liquid states, which can be 
used in numerics and experiments as indirect evidences 
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of these spin liquids. Our results can be used in further 
theoretical studies of spin liquids on square lattice, for 
example in variational Monte Carlo calculations of pro¬ 
jected bosonic spin liquid wave functions. 
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Appendix A: Solution of the algebraic PSG 

In the following we will solve the algebraic PSGs by 
using the algebraic constraints on the generators of the 
square lattice space group. 

For reasons discussed in the introduction, the low en¬ 
ergy gauge group after Higgs condensation is Z 2 , which 
is called the invariant gauge group (IGG) of the PSG. 

The two elements of the IGG are identity operator and 
the generator 

bja t bj a . (Al) 

For furture convenience, we introduce two difference 
operators A x and A v defined as A X f(x,y) = f(x,y) — 
f{x - 1, y) and A y f{x, y ) = f(x, y ) - f(x, y- 1). 

The PSG is solved under certain gauge, therefore it is 
important to consider how PSG elements change if we do 
a gauge transformation to the ansatz. 

When a gauge transformation G : bj a -A e % ^’bj a is 
applied to the ansatz, the ansatz will be invariant un¬ 
der GGxXG~ l = GGxXG~ 1 X~ 1 X, therefore Gx is 
now converted to GGxXG~ 1 X~ 1 and the phase func¬ 
tions changes according to: 

<t>x {0 -A <j>a{ r) + </>x(r) - <M^ - 1 (r)]. (A2) 

With this recipe we can use the gauge freedom to do 
the gauge fixing procedure on a spanning tree (an open 
boundary condition is assumed), thus the gauge is fixed 
to be 

<t>T! (x, y) = 0 , 0t 2 ( 0 , y) = 0. (A3) 

The commutative relation T^ 1 T 2 TiT 2 1 = 1, when 
translated into the PSG language, is 

T- [ G^G T2 T 2 G t ^T.^ [ G^ g IGG, (A4) 

therefore A x (j>T 2 (x, y) = Pi^ (in this and next section, all 
the equations are true modulo 2 - 7 t ). Here the number pi 
is either 0 or 1 due to the Z 2 gauge structure. Integers 
Pi,i = 2,3, • • • 9 appeared later are also Z 2 integers. 


After gauge fixing procedure, we are still left with three 
gauge freedom which will not change </>Ti, 4>t 2 u P to IGG 
elements but may affect other PSG elements. 


G ! : 

y) = const., 

(A5) 

G 2 : 

h (x,y) = nx, 

(A 6 ) 

G 3 : 

fa (x,y) = tt y. 

(A7) 


We are left with commutative relations, a 1 TicrT 1 1 = 
cr~ 1 T 2 aT 2 = a 2 = a~ 1 C 4 aC 4 = C\ = C , 4 “ 1 T 1 C' 4 T 2 = 1, 
and four commutative relations concerning time reversal 
symmetry, T^faTT^ 1 = T^'T 2 TT 2 l = T^aTcj- 1 = 

t- [ c 4 tc^' = 1 . 

Using the condition a~ l T\aT^ 1 = 1, and a~ 1 T 2 aT 2 = 
1 , we can obtain: 

A x ( h(x,y) = p 2 7 T, (A 8 ) 

A v (j> a {x,y) = p 3 n, (A9) 

whose solution is 4> a {x, y) = p 2 irx + p^ny + </v(0,0). 

The constraint from cr 2 = 1 is 2^(0, 0) = p 4 n. Due 
the Z 2 gauge structure, an overall phase tt has no conse¬ 
quence and hence we can fix <p a (0, 0 ) to be p 4 ir/2. 
Therefore we have 

<j> a {x,y) = p 2 nx +p 3 ny +p 4 n/ 2 . (A10) 

The gauge condition Gi, G 2 and G 3 have no effect on 
it. 

Consider the fourfold rotation G 4 . 

From G 4 T 1 TiG 4 T 2 = 1 and C^T^CrXj -1 = 1 we have 

Xxfici = Pi^y+P5^i (All) 

A y 4>c 4 = pinx + p 6 TT, (A12) 

the solution of which is (^c 4 (x, y) = pinxy+p 5 TTX+peTry+ 

0c 4 (O,O). 

From C 4 aC 4 a 1 = 1 we have 

P2+P3+P5+P6 = 0, 2 <j) Ci = P7K. (A13) 

The gauge condition G 2 can tune p$ to be zero. And 
the condition C 4 = 1 gives no new constraint. 


Therefore we have 

</>c 4 (x, y) = pmxy + p 6 iry + prnl 2. (A14) 

As for time-reversal symmetry, the condition 

T-'GTT-r 1 = 1 and r" i T 2 TT 2 _1 = 1 yield 

A x (j)T = P8K, (A15) 

Aj y <f>T = Ps'x- (A16) 

Therefore we have 

(t>r{x,y) =psTrx+p 9 TTy + (j)r{0,0) (A17) 


The condition T 1 a'Tcr 1 = 1 produces no constraint. 
The condition T~ 1 C 4 TC^ 1 = 1 requires that p 3 = pg ■ 
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Note that the gauge condition G-\ (bj a ) = e l(l, 1 bj a acts 
nontrivially in respect of the time reversal symmetry 
since time reversal operator will change a number to its 
complex conjugate. After gauge transformation G i, the 


phase function changes as 

(/> r (r) -> 0r(r) + 20i (A18) 

Therefore we may use the gauge condition G\ to fix 
</> 7 "( 0 , 0 ) to be —p$n/ 2 , so we finally obtain 

0r(r) =p s {x + y- 1/2 )tt. (A19) 

Finally the algebraic solutions of PSG are: 

= 0, (A20) 

<t>T 2 = Pikx, (A21) 

< i>a = P 2 KX + p 3 iry + (A22) 

0c 4 = PiTrxy+ (p 2 +p 3 )ny + 1 ~y, (A23) 

<t>T = P8(x + y- 1/2)tt. (A24) 


There are two choices for each pi(i = 1,2, 3,4, 7, 8 ), so 
total amounts of PSG solutions are 2 6 = 64. 


Appendix B: Physical realizations of PSG on square 
lattice 

PSG solved by the commutative relations of group gen¬ 
erators is called the algebraic PSG since it relies only on 
the structure of the symmetry group and may not be 
realized by a mean-field ansatz. 

Generally speaking, realization of a particular kind of 
ansatz on square lattice will impose further constraints 
on algebraic PSG, therefore it is necessary to analyze the 
constraints imposed on PSG when demanding that an 
arbitrary bond Uij (Uij represents A ZJ or Bij) is nonvan¬ 
ishing. 

We shall only consider bonds that start from the origi¬ 
nal point ( x , y) = ( 0 , 0 ) since other bonds can be obtained 
by translation. 


1 * + 0 , 0 ) — 


The bond A( 0)0 )_>.( x , x ) becomes A( X;X )_q 0i0) under sym¬ 
metry operation T^TfGf. Since A t j is antisymmetric, 
i.e. Aij = —Aji, we have: 


[GT 2 T 2 ] a: [GT 1 T , l] :r [Gc4C , 4] 2 A(0.0)-).(x,x) = -A(o,o)^(iz,:e), 

(Bl) 

thus we have 

2(t> Ci (0,0) + </>c 4 (~x, x) + 0c 4 (~x, -x) 


X 

+ (0, -x + i) + (j) T2 (x,i)) = n, 


(B2) 


therefore piitx 2 = n. 

A(o,o)->(a;,x) is invariant under the combined operation 
Gju, so we have: 


<^C 4 (0,0) + <pcA x -> x ) + 0) + <j> a {x,-x) 

= Pinx + (P4 + P 7 )tt = 0 . 


(B3) 


Therefore when x is odd, the condition (B2) and (B3) 
are satisfied only when p\ = p± + pj = 1 , and when x is 
even, the two conditions cannot be satisfied simultane¬ 
ously and hence A( 0 i 0 )_s.f x vanishes. 

3(0,0) 
have 


>(x,c 

>(x,x) becomes -^^ 0 , 0 )- 


(x,x) under 7?Tf G|, so we 


</>c 4 {-x,x) + (j)Ci(—x, -x) - 2 (j) Ci (0,0) 

X 

+ '52(<t>T 2 (0,~x + i) -<j>r 2 (x,i)) (B4) 

i —1 

= Pi ttx 2 = Arg (B* { 0 fi) ^ {xx) /B im ^ {X ' X) ). 

B( o,o)->(x,x) is invariant under G 4 CT, so we have 


— 0C 4 (O, 0) + (/>C 4 (x, x)—<ficr(0, 0)+(/><r(x, ~x) = pilTX 2 =0. 

(B5) 


When x is even, the condition (B4) and JB5| demand 
that f?(o,o)->(x,x) is r eal. When x is odd, ( |B5[ ) demands 
that pi = 1 and (B4) demands that £?( q 2 o)-^(x,x) is real. 
The results are summarized in Table Uvl 


-^(0,0) — >(x,x) 

-B(0,0)—>-(x,x) 

x (mod 2) 

vanish 

Im(B)=0 

0 

pi =P 4 +PT = 1 

pi = 0, Im(B)=0 

1 


TABLE IV: Constraint imposed on PSG by nonvanishing 
bond +o,o)—>(x,x) 


2 - W( 0 , 0 )->(j:, 0 ) 


The bond A( 0 ,o)->(a;,o) is invariant under a , so 

</> CT (0,0) + 0 a {x, 0) = p 2 nx +p 4 ir = 0. (B 6 ) 

This bond becomes its inverse under TfGf, thus we 


have 


2<j> Ci (0,0) + <t> Ci (0, x) + (j) Ci {-x, 0) 
= (P2 +P3)nx = 7r. 


(B7) 


When x is even, the condition (B71 cannot be satisfied. 
When x is odd, the cond itio n (B 6 | demands that p 2 + 
P 4 = 0 and the condition (B7) demands that p 2 +P 3 = 1. 


B(o t o)- 4 (x,o) is invariant under a, so 

- <^( 0 , 0 ) + <j> a {x, 0 ) = p 2 ttx = 0 . 


(B 8 ) 


This bond becomes its conjugate under TfGf, there¬ 
fore we have 


- 2</>c 4 (0,0) + 4> Ci (0, x) + (j) Ci ( -x , 0) 
= Arg (-B( 0i0 )_,.( a ; i0 ) / 3 ( 0 , 0 )->.(a:, 0 )) ■ 


(B9) 
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When x is even, we have Arg(F?*/ B)=0. When x is 
odd, we know from ( |B 8 [ ) that p 2 = 0. And the condition 
(B9| demands that Arg( B* / B) =p 3 7 r . 

The results are summarized in Table m 


A(o,o)—v(x,o) 

B( o,o)->(x,o) 

x (mod 2) 

vanish 

Im(B)=0 

0 

P 2 +P 3 = 1 ,P 4 =P 2 

P2 =P 3 = 0, Im(B)=0 
or P2 = 0 ,P 3 = 1, Re(B)=0 

1 


A(0,0)-*-(z,y) 

B( 0,0)->(x,y) 

x (mod 2) 

y (mod 2) 

vanish 

Im(B)=0 

0 

0 

P 2 T P 3 — 1 

P2 +P 3 = 0, Im(B)=0 
or vo. + Vo = 1. Re(13)=0 

0 

1 

1 

0 

Pi = 1 

Pi = 0, Im(13)=0 
or pi = 1 Re(B)=0 

1 

1 


TABLE VI: Constraint imposed on PSG by nonvanishing 

bond tT(0,0)— l(x,y) 


TABLE V: Constraint imposed on PSG by nonvanishing 
bond M(0,0)-Ka:,0) 


The zero energy eigenvector of the mean-field Hamil¬ 
tonian (22) at Q = ( 7 t/ 2 ,— 7 r/ 2 ) point is 


3. U(o,0)—>(x,y)j x 7^ V 


&1 




(Cl) 


There is only one constraint: 
[G ra r 2 ]"[G ri T 1 ] !E [G 04 C74] 2 A (0i0) _ KxiI , ) = -A m ^ XiV) , 
which leads to 


20c 4 (0,0) + <\>c 4 (~y, x ) + <t> Ci {-x, -y) 
y 

+ '^2(fa 2 (0,~y + i)+fa 2 (x,i)) (BIO) 

i =1 

= PlTTXy + (p 2 + Pz)Tr{x -y)= 7 T. 


When x and y are both even, the condition (BIO) can¬ 
not be satisfied, therefore A( 0 ,ow{a,y) vanishes. When 
x + y is odd, the condition (BIO) is satisfied when 
P 2 + P 3 = 1. When x and y are both odd, the condi¬ 
tion (BIO) demands that p\ = 1. 

As for B( 0 ,o)-Kx,y), we have 


- 2</>c 4 (0,0) + (\> Ci {-y, x) + (j) C4 {-x, -y) 
y 

+ ^{fa* (°’ ~y +*) - fa 2 Og *)) 

2—1 

= piTTxy + (j>2 + P 3 )tt(x - y) 

= Arg(B* 0)0 )_ > ( a!) j,)/S(o i o)_>(a,K))- 


(BH) 


When x and y are both even, the condition (B10) 
demands tha t -Bfo .Qizi's.id is real. When x + y is odd, 
the condition B10 demands that Arg (B* / B)=( p 2 +p 3 )tt- 
When x and y are both odd, the condition ( B10| ) demands 
that Arg(B*/B)=p 1 n. 

The results are summarized in Table Ivll 


When boson condense, a ground state expectation 
value at Q is obtained: ('F(Q)) = z’Fi. 

The zero energy eigenvector at — Q point is 

) ■ (C2) 

\/Bf +A^-bJ 

And the condensation value is (\H(—Q)) = 2 . 

Therefore we can represent the condensation on lattice 
site r as: 




_ — iAi _ 

VBf+Al-B ! 
e -iQ r 



(C3) 

The 2x2 matrix is propotional to a SU(2) matrix and 
its effect on the order parameter (S(r)) is just a SO(3) 
rotation. 

The magnetic order obtained from the spinor x there¬ 
fore can be computed using 


(S(r)) = (1/2 )x^ax 

l _|_ ty+y ^ i _ ^ (C4) 

2 S 2 nv ' 1 

If we introduce the notion s = , 7 lA l -, the vector 

y B \+A\-B 4 ’ 

n u and n v can be conveniently represented as 

n u = (0, 2*s, —1 — s 2 ), (C5) 

n v = (0, —2 is, —1 — s 2 ). (C 6 ) 


Appendix C: Derivation of the magnetic order from 
the condensation of Schwinger boson in mean-field 
level 


And when s / 0, we have n u ■ n v /(\n u \ ■ |n u |) = 
—1> —1. Therefore the magnetic order obtained 
from the Z 2 [ 0 , 0 ] spin liquid is a non-collinear order. 


1. Canted Neel order from ^[0,0] state 


2. Neel order from ^ 2 [0,7r] state 


For Z 2 [0, 0] state, the spinon condensation happens at 
±Q, Q = (tv/2, — 7 t/ 2 ) when p = 2B 2 — 2^/A\ + Bf. 


For Z 2 [0,7t\ state, the spinon condensation happens at 
±Q, Q = (n/2, 7 t/2 ) when p = —2Ai — 2B 2 . 
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The zero energy eigenvector at the Q = ( 7 t/ 2 , 7 t/ 2 ) 
point is 

= Q' (C7) 

When boson condense, the condensation value at Q is 
WQ)> =z*i. 

The zero energy eigenvector at — Q point is 


To = 


1 


(C 8 ) 


And the condensation value is (\I/(—Q)) = wf 2 - 
Therefore we can represent the condensation on lattice 
site r as: 


x = ( \ _ ( *, w 




vlW\ —iz 


e iQr 

o-iQ r 


(C9) 


where the 2x2 matrix is propotional to a SU(2) matrix 
and its effect on the order parameter (S(r)) is an SO(3) 
rotation. 

The magnetic order obtained from the spinor x is 
therefore 


(CIO) 


(S(r)) = (1/2 )x^ax 
= cos(2Q ■ r)m = (—) x+!/ m, 

which is the Neel order. 

3. Magnetic order from 7r-flux states 


For 7 T flux state with only nearest bond A i, the spinon 
condensation occurs at ±Q, Q = (tt/2,tt/2) when y = 
-2 A-!. 

We devide the system into two sets of sublattice which 
are distinguished by (—l) y which are labeled by u and v 
separately. 

At Q, the Hamiltonian has two linearly independent 
eigenvectors: 


i )T 


(Cll) 


The condensate at Q is the linear combination of the 
two eigenvectors 


(T(Q)) =z 1 'f 1 +z 2 ^ 2 - 
At —Q point, the two eigenvectors are: 




(C12) 


( 013 ) 


The condensation at Q is the linear combination of 
H /3 and '£'4 

(T(-Q)) = Wi ^3 +W 2 ^4- (C14) 

In real space, we have 

({Kvt)\ 

(bvrt) = e iQr^ i + ^ 2 t 2 ] + e” iQ ' r [wi'F 3 + u 2^4]- 
WurJ./ 

(CIS) 

Therefore the condensate on lattice r can be repre¬ 
sented as: 


T = ( (b U rt) 

" [ (Krl) y 


~ 1 ’ - 75 ^ 1 - 75^2 

75* 75* 2 


e'Q r 

o-*Q-r 


(C16) 


” {(bvri), 

( * 1 , 


72 


Wl 


72 

1 


W 2 


W w i’ 75 z i~75 z 2 


e'Q r 

r'Qr 


(C17) 


In order to have uniform magnitude for the ordered 
moments (l/2)x' u v x u , v , we have to impose the constraint 


x' u x u = x],x v = const. 


From (C16l and (C17) we find that there are four dis¬ 


tinct spinors in a plaquette which are named as Z\ if x 
is even and y is even, Z 2 if x is odd and y is even, Z 3 if 
x is odd and y is odd and Z 4 if x is even and y is odd. 
Only two of them are linearly independent. 

Therefore we choose two linearly independent spinors 
as: 


tp = Z 1 = ^ 
i x = Z 2 = i 


^ - 75 Wl - 75 u ’ 2 

sk T >k 1 ^ 

- 75*1 - 75+7 

( -i + 75 w i + 75 w 2 


-w? + 


75 1 


1 * 
VS* 2 , 


Therefore Z 3 and Z 4 can be represented as: 


Z» = 


z 2 + 71^1 + 
1 


75 


w 2 


\~ w 2 + Tf-l + 71 + 
= -\f2ioyty* + ix, 


Za = i 


* 2 “ 75^1 - 75 W2 

V -^2 - ~ 3 = 2*1 

= V^icTyX* + itp. 


(C18) 

(719) 

(C20) 

(721) 

(C22) 

(C23) 

(C24) 
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The constraints that condensate on every site is uni¬ 
form require that 

W 2 = l*x| 2 = I - V2icr v ip* + ix | 2 = | V2ia y x* + iip\ 2 - 

(C25) 

The constraints can be simplified if we define two 
spinors 


<t>+ = V% + + VyX*)> (C26) 

<f>~ = V2-V2(*p-a v X*)- (C27) 

Therefore we have: 

Zx = ^(^2-^0+ + V 2 + v^_), (C28) 

Z 2 = 2 ^-*a y ( ^2-V2tf + + ^ 2 + v / 2<^)(C29) 
Z3 = 2 + - V 2 - V / 20-)(C3O) 

Z4 = ^(V 2 + ^+- \/2-V2<t>-). (C31) 

And the constraints now become 


|0 +| 2 = |<M 2 , Re(40-)=O. (C32) 


Define vector order parameter on the four sites as n$ = 
ZjaZ,. (i = 1,2,3,4). It is easy to verify that Yh™i = 0. 

i 

Therefore we can parameterize vector order parameter 
on every site as: 

n = o) + {~) y n(o,n) + (~) x+y n(n,n)- (C33) 

The three vectors, n^g), ^(o,tt) and corresponds 

to mi, m 2 and m 3 in main text respectively. 

The three order parameters nr^, 0 ), n( o,tt) and 
can be represented with two spinors (p+ and cj>_ 


^(tt,o) = ((™i + « 4 ) - («2 + n z ))/A 
= \{n^ + +« 0 _), 

= (Oh + n 2 ) - (n 3 + n 4 ))/4 

= ^ (</>+#- + <?<£+)> 

^(tt,tt) = ((«i + ^ 3 ) - (n 2 + n 4 ))/4 

= -g-(-«■*>+ + »v_), 


(C34) 


where we have defined n^ ± = </>j_ B(j)±. 

The manifold of the order parameter can be described 
as follows. Given two spinors <f>+ and </>_ with the 
same magnitude, we can perform SU(2) transforma¬ 
tions to them separately without altering the condition 
\4 >+\ 2 = |</>-| 2 - And then we can tune the relative U(l) 
phase between the two spinors to meet the condition 
R e{(j^ + (j>-) = 0. 


In principle, the transformation rule of the spinor fields 
< t>+ and </>_ under space group operation can be worked 
out through the corresponding transformation rule of 
bosonic field under PSG. 


It is easy to deduce from constraints (C25) that the 
intersection angle between magnetic moments on neigh¬ 
boring sites should always be greater than 90°. And note 
that the Neel state is not forbidden by this condition and 
therefore is one of possible magnetic ordered states that 
can be obtained from the 7r-flux state through boson con¬ 
densation. 


Appendix D: Derivation of the continuum-field 
theory for the transition from Z 2 spin liquid states 
to Neel ordered states 

In this Appendix we derive the continuum-field theory 
for the two Z 2 spin liquid states from the microscopic 
Hamiltonian in the long wavelength limit. 

We divide sublattice u and v with different parity 
(—) x+y : u sublattice is parity-even and v is parity-odd. 

Following Read and Sachdev’s prescriptiorP^, we repre¬ 
sent the boson operator in terms of the two slow-varying 
fields 


h — ib irlc^^ 1 

u u,r,a. 'ru,<y.\ L •> 

b v , r , a = -cr^* i/3 (r)e iQr , 


(Dl) 


where Q represents the k point where condensation of 
bosons happens. 

The bond operator Ai, B\, B 2 therefore can be repre¬ 
sented as: Nearest bond A 1 : 

4fr,r' — 2 bu,r, 

= + Ar • 9 r + (Ar 2 9r)2 + . ..]<«• 

(D2) 


Nearest bond B\: 


B r ,r' — 2 ^u,r,abv,r',a 


= -U Q Ar < 0 r u , a \ 1 + Ar . a r + 


• ••RV 

(D3) 


Nex-nearest bond B 2 : 


1 6 t b , 

2 u u,r,a u u,r,ai 

= \e iQ Ar r u , a [ 1 + Ar . a r + (Ar 9 5r)2 + . • .]Wa, 


(D4) 


2 
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and 


b / 

2 u v,r,a u v,r' ,a 


= \e iQ Ar r v ,a[ 1 + Ar . 9 r + (Ar 2 9r)2 + .. .]i, ai 

(D5) 


where Ar = r' — r. 


1. Z 2 [ 0 , 0 ] state 


£= d 2 r{- - i on -i 2 d T z* a -d T z a 

J (—n + IA\ + 2tf 2 )a 


+ (Ai - 2 B 2 + 


4|Bi 


~H + 2 Ai + 2B2 
, — 2 /i - 4A 1 + 4B 2 


)3 r z* • <9 r z Q 


(—/r + 2Ai + 2 I? 2 )a' 


\Bi?)zlz 0 




(D9) 


The minima are located at ±(7t/ 2, —7t/2). Therefore 
we choose Q = ( 7 r/ 2 , — 7 r/ 2 ) in (D1). 

Inserting (|D2|)-(D4) into the mean-field Hamiltonian 


( 22 ), we obtain: 


C = 


J ^{C,a ~ ^ ~ MC.aVV 


^v,cA^v,a) 2,Aiip u ^ a i!) v a T • d r 'i/j v a T c.c. 


~ 2B l i ° y a fji l l, 0 .' l l J *v,p + y B *i^fid r ipl ta ■ d r + c.c. 

+ 2B 2 {’lpl ta ^ u ,a + VC.aV’w.a) ~ Cl 2 B 2 (d T lj)* ua • <9 r V>u,a 

+ 5 r l/>* ia • d r V!„, a )}. 


(D 6 ) 


Note that terms with odd spatial derivatives vanish 
due to the geometry of the square lattice. 

Introduce two fields: 

Za (l^u,a T VV,a)/2> Pu,a '*Pv,a')/^- (D7) 

The Lagrangian now becomes 


Note that the Higgs term B\-^z a ia v a pZp + c.c. plays 
the role of the vector spin chirality order parameter for 
reasons discussed in the text. 

The Lagrangian can be further simplified after in¬ 
troducing two fields ui a , a = 1,2 (here we assume 
Im(I3 1 )>0), 


w 1 


w 2 


zt + z i + i{z } - zp 
2 

4 + z} + i(z t - Zj.) 

2 


(DIO) 


and has the form 


A = / d 2 r{- - , , 2 9 t u* • d r w Q 

J (— f-i 2Ai H- aB2)cl z 


+ (Ai - 2 H 2 + 


4|Hi| 2 


— /U + 2 Ai + 2B2 
, -2 n - 4Ai + 4I3 2 


)d r w* • a r w 0 


(—/r T 2Ai + 2B 2 )a 


(—/r + 2 Ai + 2 I? 2 )a 2 

Im(Hi)w* ^-w a }. 
z ar 


|Si| 2 )o;> a 


(Dll) 


_ f d 2 r ro ^ d * d 

— / 9 7 ^Ta: ”1“ ^TTo; j -2a: 

J a z ar ar 

+ ( — 2/i — 4Hi + 4B 2 )z*z a + (—2/r -f- 4Hi + 4 .B 2 ) 7 t* 7 t 0 , 
+ a 2 (Hi - 2 B 2 )d r z* a ■ d r z a 

- a 2 B\i(j y a pdrZ* a ■ drTTp + c.c. + ABlz* a ia v a p'K*p + c.c.}, 

(D 8 ) 


Therefore the Lagrangian will flow to a new fixed point 
where space and time scales differently. The dispersion 
relation is now ui oc k 2 with dynamical critical component 
z = 2 . 

The term with quadratic powers of time derivatives is 
irrelevant and hence can be dropped. 


where terms involving it field and spatial derivatives are 
omitted since they will generate terms with fourth or 
higher power of spatial derivatives of z field after inte¬ 
grating out 7 T field. 

Fields 7r a have a large mass gap — 2/z + AAi + 4B 2 and 
can be safely integrated out. The low energy effective 
Lagrangian after integration is : 


2. Z 2 [0, 7r] state 


In this case the minima of the spinon dispersion are 
located at ±Q where Q = {tt/2, 7 t/ 2 ). 


Inserting (D2)-(D4) into the mean field Hamiltonian, 
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we therefore obtain: 


r d 2 v d d 

£ = I ~^2~{' l l ) u,a~^ , J’u,a ~ 1pv,a V’t'.a — A t (V , M,a’ i / , M,a 

a 2 

+ VC.aVV.a) - 2+ —Axd^u^oc ■ d r ^* a + C.C. 

+ y^< / 3 C,J9 2 -5 2 ]^ + c.c. 

- 2- B 2«, Q V , «,a + C.aV^a) + a 2 £ 2 (3rC,a ■ d r 1p u ,a 

+ d r 1PI a ■ dr 1p v ,a)}. 

(D12) 


Introduce two fields 2 and 7 r defined in (D7), the La- 
grangian becomes: 



FIG. 13: Gauge choice on the dual lattice (solid black lines). 
The vison hopping Jij on the thick bonds are negative. The 
blue dashed line represents the original lattice. Note that all 
the point group symmetry operations are defined with respect 
to the original lattice site O. 


r - f d2r so * d _1 n * d 

C ~ J 0? {2Za dT 7Ta+2na dT* a 
+ (—2/i - 4Ai - 4 B 2 )z*z a + (- 2 /X + 4 A x - 4R 2 )<7r Q 
+ o?{A\ + 2 B 2 )d r z * • d T z a 
- a 2 Blia y al) z* a [dl - d 2 ]w* p + c.c.}. 

(D13) 

After integrating out 7 r field, we obtain a low energy 
effective Lagrangian: 


Following Ref. 031 the dynamics of vison can be de¬ 
scribed by a fully frustrated transverse-field Ising theory 
on the dual lattice. 

H = Y, J n T i T j-Y, K iTi--- > ( E1 ) 

ij i 

where the product of bonds around each elementary pla- 
quette is 


C = / d 2 r{ 


(—/r + 2 Ai — 2B2)' 


d T z* ■ d T z a 




+ (Ai + 2 B 2 )d r z * • d r z a + 


(-2jt - 4Ar - 4R 2 ) 


(D14) 


Here the Higgs terms consist of B\ have cubic powers 
of time and spatial derivatives or higher and hence is 
irrelevant by naive power counting. 

After rescaling and restore the compact gauge field, we 
obtain an effective field theory cosists of a massive boson 
z a coupled to a compact U(l) gauge field, 


C = / d 2 r[|£>. 


■ z\ 2 + c 2 \D r 


■m 2 \z\ 2 ] 


(D15) 


The critical point of this theory is p = — 2Ai—2R 2 consis¬ 
tent with the mean field solution and the spinon velocity 
c is proportional to ^Ai{Ai + 2 B 2 ). 


Appendix E: Calculation of vison PSG 

In this section we will derive the vison PSG for square 
latticd^®. 

The dynamics of vison can be described by an odd Ising 
gauge theory in the dimer limit. After a duality trans¬ 
formation, the odd Ising gauge theory is transformed to 
a transverse field Ising model on the dual lattice. 


n s S n ( J b) = -!• (E2) 

plaquette 

This Hamiltonian is invariant under Z 2 gauge trans¬ 
formation 


Ti t I]/ Tj . K, t 7/, K t . t fj r Tjj Jjj . (F3) 

where r)i = ± 1 . 

Due to the Z 2 gauge structure, we adopt a specific 
gauge choice (FIG. 13 1 for J,.j and calculate the vison 


PSG under this gauge. 

In the following we shall ignore the kinetic term K)rf 
and adopt a soft-spin formulation where the vison field 
rf’s take real values. Considering nearest and fourth 
nearest neighbor interaction of the vison fields as consid¬ 
ered in Ref. 03J there are eight (rather than 4 because two 
sets of sublattice are not equivalent) inequivalent mini¬ 
mal points in the Brillouin zone (see Fig. 14). 


Qi = ( 0 , 0 ), Q 2 = ( 0 ,7r), Q 3 = (-, —), Q 4 = (-, -), 

7r —7r 7r 7r 

Qv = ( 7 T, 7 r), Q 2 > = (tt, 0 ),Q 3 / = (-, — ),Q 4 ' = 

(E4) 

Expand the vison field with slow varying modes 
(j> a (a = 1, 2, 3,4, V , 2', 3', 4') at these eight momenta: 


= £< 


JQa.-r 


(E5) 
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FIG. 14: Brillouin zone of the square lattice, and the eight 
momenta, labeled as 1, 2, 3,4,1', 2', 3', 4' and indicated by red 
squares in the figure, for low energy vison modes. 


Under symmetry operation, the vison fields transform 
as 


A : <v 

+> 

n x+1 T : + i^ 


(E 6 ) 

T 2 ■ T X: y 

-+ 

(-)*<„+!, 


(E7) 

P • T Z 

r x • I x ,y 

-+ 

(_\x_z 
\ ) ‘x,l—yi 


(E 8 ) 

P • T Z 

r y • ‘x,y 

-+ 

T Z 

l—x,yi 


(E9) 

P ■ T Z 

-+ 

(1 - i)i x+ y + (1 + i)i- x -y _ z 

.(E10) 

r xy • > x ,y 

2 

‘y,x 

We choose 

base 

functions 

as 

(<A, 02, 4>3, 04, 01', 02', 03', 04') T , 

the vison 

modes 


transform under smmetry operation as 

8 

0n t ^ ^ [^i^(+)] n ,m0m- (Ell) 

m—1 


Therefore we can directly write out the transformation 
matrices 


<V!Ti) = 


(0 

0 

0 

0 

0 

1 

0 

v° 


0 

0 

0 

0 

1 

0 

0 

0 


0 0 0 
0 0-1 
0 i 0 
-i 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


-1 0 o\ 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 -i 
0 i 0 / 


(E12) 


/o 0 0 0 0 1 0 o\ 

0000-100 0 

0 0 0 i 0 0 0 0 

0 0 * 0 0 0 0 0 

0 -1 00000 0 
1 0 0 0 0 0 0 0 

0000000 -i 

\0 0 00 0 0 -i0/ 


o^{P x ) 


/o 0 0 0 0 1 0 o\ 

00 0 0-1000 
0 0 0 0 0 0 *0 

0 0 0 0 0 0 0 * 

0-10 0 0000 
1 0 0 0 0 0 0 0 

00 -i 0 0000 

\0 0 0 -i 0 0 0 0 / 


(E14) 


O^Py) 


/lOOO 0 0 0 o\ 
0100 0 0 00 
0 0 0 -i 0 0 00 

0 0 * 0 0 0 0 0 
0000 -1 0 00 
0000 0-100 
0 0 0 0 0 0 0 * 
\0 0 0 0 0 0 -i 0 / 


(E15) 


0<f,{P X y) — 


/ 0 0 

0 0 

0 1 + i 


1 + i 0 
0 0 
0 0 
0 1 -* 
\l-i 0 


0 

1 — * 

0 

0 

0 

1 + * 
0 
0 


1 — * 
0 
0 
0 

1 + i 
0 
0 
0 


o o o i + *\ 

0 0 1 +i 0 

0 1 — i 0 0 

1 — * 0 0 0 
0 0 0 1 — i 

0 0 1 — * 0 
0 1 +i 0 0 

1 +i 0 0 0 / 

(E16) 


In the following we may identify the operation P x 
with er and P xy P x with C 4 , therefore 


o*(t 2 ) 


(E13) 
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0^(C 4 ) = 04P xy )0t{P x ) = - 


( 0 0 0 1 -i 0 0 0 

0 0 1 -i 0 0 0 1 + i 

1 -i 0 0 0 —1 — i 0 0 

0 - 1+i 000 1 + i 0 

0 0 0 -1 -i 0 0 0 

0 0 -1 -i 000 -1+i 

1+i 0 0 0 -1 +i 0 0 

V 0 -1-i 0 0 0 1 — i 0 


1 + i \ 
0 
0 
0 

-1 + i 
0 
0 

0 ) 


(E17) 


r 


The vison PSGs now can be easily obtained: 

0 4 ,(T2)- 1 04,(T 1 )04 > (T2)04,(T 1 )- 1 = - 1 , 
O^a^O^O^O^T,)- 1 = -1, 
0^(a)- 1 0 < ^(r 2 )0^(a)0^(T 2 ) = 1, 

OAa) 2 = 1 , 

’ (Eli 

O^a)- 1 0^)0^)0404) = 1 , 

04C 4 ) 4 = - 1 , 

04C4)~ 1 04T 1 )04C 4 )04T 2 ) = - 1 , 
04C4)- 1 04t 2 )04c 4 )04t 1 )- 1 = 1 . 

Appendix F: Solutions of Abrikosov fermion PSG 


The fermionic PSGs are characterized by four num¬ 
bers: r) X px , r] xp y , T] xy , 7j tl and numbers from the commuta- 


tive relations of SU(2) matrices as defined in TABLE VII 


Algebraic Identities 

Z 2 numbers from commutative relations 

cr 2 = 1 

[gp v ] 2 = w 

cr~ 1 C4aC4 = 1 

[9py]~ 1 [gpx y 9p y }[gp y ][gpx y gpy\ = +rc 4 

Cf = 1 

[gp^gpyV = VC 4 

c4t~ 1 c 4 t = 1 

lgPx v gPy]~ 1 gT 1 9Pxy9Py97' = VO 4 T 

a~ 1 T~ 1 aT=l 

gp^T 1 9 p y9T = +tT 


TABLE VII: Z 2 numbers from the commutative relations 
between SU(2) matrices. Here rji are numbers that can be 
±1, and g-i are SU(2) matrices as defined in Eq. (FI)-(F6). 


The Abrikosov fermion PSG on square lattice are stud¬ 
ied in great detail in Ref. [27] In this Appendix, we briefly 
summarize the algebraic solutions of fermionic PSGs and 
introduce some notations for convenience. 

The solutions to the fermionic PSGs are as follows: 


GrJ*) = vlYo, (FI) 

Gr y (i) = r 0 , (F2) 

G P*(*) = VxpxVll'pygP^ (F3) 

G Py (i) = V l x x P yV l x v p X 9P y , (F4) 

Gp. y (i) = vlf9P xy , (F5) 

G r {i) = V t x+ly 9r ■ (F6) 

In Wen’s notation, T x , y stands for translation symme¬ 
try, and three parity operations P x , P y and P xy acts as 
follows: 

Px '• {ixpy) t ( i Xl iy)i (F7) 

Py ■ (tjEjty) ^ ( ixt ’iy)i (F 8 ) 

Pxy ■ (+1 by) ^ (+j : A;) ■ (F9) 


Appendix G: Physical realization of Abrikosov 
fermion spin liquid states 

In this Appendix we will analyze the physical realiza¬ 
tion of the fermionic PSGs as shown in TABLE Mil As 
discussed in Ref. 123 the mean-field ansatz U l:j should 
take the form of 


U iJ = ip ij Wi j , (Gl) 

(here Pij is a non-negative real number, and If), £ 
SU( 2)) in order to describe a spin-rotational symmetric 
fermionic spin liquid state, and the time-reversal symme¬ 
try T is implemented as 


stands for reflection along y axis, reflection along x axis 
and respectively. 

Note that C 4 = P xy P y , therefore we can obtain the 
PSG for C 4 from the product of P xy and P y . And the a 
operation is just the operation P y in Wen’s notation. 


7” : Uij —> —Uij. (G2) 


The symmetry operations are defined in Eq. (F7)-(F9). 
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1. Realization of the ZzBtLtLt 1 ^ state 


Therefore we have 


PSG elements of this state are: 


G x (i) = (-)^ro-, 

(G3) 

G y (i) = r 0 ; 

(G4) 

Gpji) = (-)*”+** in; 

(G5) 

Gp y (i) = (-) il+i nr i; 

(G 6 ) 

Gp xy {i) = (-) ixi nri; 

(G7) 

Gr(i) = (-)’ I+i dr 3 . 

(G 8 ) 


There are two sites in a unit cell of 7 r-flux fermionic spin 
liquid states distinguished by parity (—) x = ± 1 , which 
are labeled by u and v respectively. 

First, let us consider the Lagrangian multiplier a l 0 . 
Since the term OqT; should be invariant under Gq-T, 
we have {oq 77 , 73 } = 0. And it is also invariant under 
Gp x P x ,Gp y Py and G Pxy P xy , so we have [a() 77 , 1 - 1 ] = 0. 
Hence we have aj ^ 0 and < 2 q ’ 3 = 0. 


U x +y — Uo,0+x-\-y — Go—x—y,0 — U x -\-y, (G16) 

where the first and the third equation are due to 
Eq. (G13), and the second is due to Eq. (G15) and the 
symmetry of the ansatz Ui :] = LH. 

Finally we have u x +y = u x -y = 0. 


c. Third neighbor bond 

The general form of the third neighbor bonds is: 

Ui,i+2x = tt2£, (G17) 

Ui,i+2y = U2y ■ (G18) 

Considering time-reversal transformation T, we imme¬ 
diately obtain {u 2x ,gp} = {u 2y ,gr} = 0- Considering 
reflection P y , we have [u 2x ,gp y \ = [u 2 y,gp v \ = 0. And 
under reflection P xy , we have u 2y = iTiu 2x (iriy. 

Therefore we have u 2x = u 2y = x t i- 


a. Nearest-neighbor bond 

With translation along x and y direction, we can write 
down the general form of the nearest bonds: 

Ui,i+x — u x , (G9) 

U it i+g = (-) ix Wy, (G10) 

where u x and u y are site-independent. 

Since u x is invariant under time-reversal operation, we 
have Gj-(i)u x G!j-(i + x) = — u x , which leads to [%, r 3 ] = 
0. 

u x is also invariant under P y : Gp y (i)u£Gp^(i + x) = 
u$, therefore {u£,ti} = 0. Hence we have 

Ux = VT3. (Gil) 

Under P xy , G Piy {i)u x Gp {i - l-i) = (-) ly u y , therefore 

u y = -r]T 3 . (G12) 

b. Next-nearest-neighbor bond 

The general form of the bonds is: 

Ui,i+x+y = ( ) x U x -\- y , (G13) 

Uiq + i-y = {-)^U X -y. (G14) 

Under P x P y , we have 

Gp x P x Gp y PyU 0 ,o+x+y[Gp x P x Gp y Py\ ] 

- tt 

where O represents the original point. 


d. Fourth neighbor bond 

The fourth neighbor is required to obtain a Z 2 spin 
liquid. The general form of the fourth neighbor is 


Ui^i-\-2x-\-y = ( ) x 'U j 2x+y-> 

(G19) 

Ui : i-\-2x—y = ( ) x 'U j 2x—yi 

(G20) 

Ui,i+x-\-2y = 'U j x-\-2y-> 

(G21) 

Ui^i-\-x—2y — 'U j x—2 y- 

(G22) 


Considering time reversal transformation, we have 
[u 2x +y,gr] = 0. 

For convenience we choose u 2 £+y to be ryro, which 
ensures this state to be a Z 2 spin liquid state. The other 
bonds can be obtained under reflection P x , P y , and P xy . 


The results are 

Ui,i+ 2 i+y = (-)* x * 7 T o, (G23) 

Ui, i+ 2x-y = -(-^Upo, (G24) 

Ui, i+x+2 y = rpo, (G25) 

Ui, i+x - 2y = -i'jTo. (G26) 

In conclusion, the symmetry allowed ansatz for the 
Z 2 Bt}_t}_t 1 t^ state are 

Ui, i+X = V t 3 , (G27) 

U^ +ij = -(-) ix r?r 3 , (G28) 

Uiq+ 2x = Uiq+ 2 y = X T h (G29) 

Ui,i +2x+ y = (-)*n 7 To, (G30) 

Ui,i +2x -y = -(-)* X * 7 A>, (G31) 

Ui,i-\-x+2y = *770, (G32) 

U iji+x - 2 y = -ryr 0 , (G33) 

aj + 0. (G34) 
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After Fourier transformation, the Hamiltonian be¬ 
comes H = where 

M = 2r/cos(k x )Ti — 2r]COs(k y )T 2 (G35) 
+{aj + 2x[cos(2fc x ) + cos(2fc y )]}r 3 (G36) 
+ 27 [sin(fc x + 2k y ) — sin(fc x — 2k y )]T± (G37) 

+ 27 [sin( 2 fc x + k y ) — sin(2fc x — fc y )]F 5 , (G38) 

where 77 , x, 7 are all real numbers, and we have used the 
Nambu spinor representation 


= (/«kt, fv kt, fL kt, (G39) 

and 

ri=r 3 ( 8 iTi, (G40) 

r 2 = T 3 (8)r 3 , (G41) 

F 3 = n 0r o , (G42) 

r 4 = r 0 (8)r 1 , (G43) 

r 5 =r 0 (g)T 3 . (G44) 


Note that k x € (0, 7 r), k y € (— 7 r, 7 r). 

When there is only nearest-neighbor bond present, 
the state will have two Dirac points at wave-vectors 
(A : x ,k y ) = (7t/2,±7t/ 2). By including the Lagrangian 
multiplier aj, the third neighbor bond and the fourth 
neighbor bond, we find that two gaps open at the two 
Dirac points. 


2. Realization of the Z 2 Br!r!r 1 d state 
PSG elements for this state are: 


G x (i) = (-) iH T 0 ; 

(G45) 

G y {i) = r 0 ; 

(G46) 

Gp x {i) = {—) ix+iy iT 2 i 

(G47) 

Gp y (i) = (-) ix+iy iT 2 i 

(G48) 

= (-) y dr 1 ; 

(G49) 

G r{i) = (-) ix+i nr 3 . 

(G50) 


First, since the term a l 0 Ti should be invariant under 
Gj-T, we have {a l 0 Tp r 3 } = 0. And it is also invari¬ 
ant under Gp x P x and Gp P xy , so we have [oqT^ti] = 
[a l 0 Ti,T 2 ] = 0. Hence the Lagrangian multipliers aj ' 2,3 = 
0 . 

The analysis of the symmetry allowed ansatz is similar 
to the previous case. A possible ansatz for the realization 
of a gapped Z 2 spin liquid state is 


Ui,i-\-x — V r 3, 

(G51) 

Ui,i+y = -(-)^r/r 3 , 

(G52) 

Ui t i+2x = X r 2 > 

(G53) 

Ui,i+2y = — XT 2 , 

(G54) 

^i,i-\-2x-\-y = ( )" x i / yTo 5 

(G55) 

Ui t i+2x-y = -(-)" x i7^0, 

(G56) 

Ui,i-\-x-\-2y — ^7^0? 

(G57) 

i,i-\-x—2y = ^7 r 0? 

(G58) 

1,2,3 n 

a o = 

(G59) 

After Fourier transformation, the Hamiltonian be- 

comes H = 4 , j c A/\E , k. where 


M = 2?ycos(fc x )ri — 277cos(/c y )r 2 
+ 2 x[cos( 2 fc x ) — cos( 2 /c. y )]r 3 
+ 27 [sin(fc x + 2k y ) — sin(fc x — 2 fc y )]F 4 
+ 27 [sin( 2 fc x + k y ) — sin( 2 /c x — fc y )]F 5 , 

(G60) 

where 77 , 7 are all real numbers, and we have used the 

Nambu spinor representation 


’hk = (f u k|) fv kt> fu~ kf fl- 14) ’ 

(G61) 

and 


Ti = t 3 <8 Ti, 

(G62) 

r 2 = r 3 0 r 3 , 

(G63) 

r 3 = t 2 0 r 0 , 

(G64) 

r 4 = T 0 O Ti, 

(G65) 

r 5 = r 0 <8 7 - 3 . 

(G 66 ) 


Note that k x € ( 0 , 7 r), k y € (— n,n). The fourth neigh¬ 
bor bond opens up gaps at two wave-vectors (k x ,k y ) = 
(7t/2,±7t/2). 

3. Realization of the t£ state 

PSG elements are: 


G x (i) = (-) iH T 0 ; 

(G67) 

Gy(i) = r 0 ; 

(G 68 ) 

Gp x (i) = (-)^ +i nri; 

(G69) 

Gp y {i) = (-) ix+i nri; 

(G70) 

Gp iy («) = ( —)* x h/ To ; 

(G71) 

Gr(i) = (_ y*+ l yi T3 . 

(G72) 


The analysis of the symmetry allowed ansatz is simi¬ 
lar to the previous case, hence we will directly show the 
results here. The next-nearest-neighbor bonds are again 
prohibited by the symmetry in this case. 
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The symmetry allowed ansatz for the ^-Bririr 0 ^ 
state are 


Ui,i-\-x — 77 T 3 , 

(G73) 

Ui,i+y = (-)*» 777-3, 

(G74) 

2x = 'U j i,i-\-2y = 'XTl 5 

(G75) 

Ui_ i i-\-2x-\-2y ~ ^2 1 

(G76) 

Ui,i-\-2x—2y — ^^~2i 

(G77) 

aj^0. 

(G78) 


After Fourier transformation, the Hamiltonian be¬ 
comes 

H = 2rjcos(k x )Y i + 2ijcos(k y )T 2 
+ {oj + 2x[cos(2fc x ) + cos(2fc y )]}r 3 (G79) 

+ 2£[cos(2k x + 2 k y ) — cos(2k x — 2 fc y )]r 4 , 

where k x £ (0, n), k y £ (7r,7r), and 

ri=T30r 1; (G80) 

r 2 =T3®T 3 , (G81) 

r 3 =r 1 ®r 0 , (G82) 

r4=T2®T 0 . (G83) 

The Lagrangian multiplier aj and the third neighbor 
bond open up gaps at the two Dirac points ( k x ,k y ) = 
(tt/2, ±71-/2). 

4. Realization of the Z 2 AtLt-T°t + state 
PSG elements are: 


G x (i) = G y (i ) = r 0 ; 

(G84) 

Gp.(i) = (—)**+*» in; 

(G85) 

G Py {i) = (—) ix+iy iTT, 

(G 86 ) 

Gp xy {i) = To; 

(G87) 

G r (i) = ir 3 - 

(G 88 ) 


First, since the term ci 1 0 ti should be invariant under 
G'fT', we have {cl 1 qTi, t 3 } = 0 . And it is also invariant un¬ 
der Gp : P x ,G Py Py and G Pxy P xy , so [a l 0 T U Ti] = 0- Hence 
we have aj 7 ^ 0 and Oq ’ 3 = 0 . 

a. Nearest-neighbor bond 

In this state, the gauge transformation G x ,G y: G Pxy 
are trivial, we can simply write down the general form of 
the nearest bond: 

— u x . (G89) 

The bond u x is invariant under time-reversal transforma¬ 
tion, therefore we have G'p{i)u x G'p{i + x)^ = —u x , which 
leads to {u x ,r 3 } = 0 . 


u x is also invariant under P y , therefore we have 
Gp y u$Gp = u x , which leads to {it£, ti} = 0 . 

In conclusion, the nearest bonds are as follows: 

Ui,i+x — u i4+y = tjt 2 . (G90) 

b. Next-nearest-neighbor bond 

The general form of the next-nearest bonds are 

“ tlx+y, (G91) 

Gi'i-\- x —y t^x—y- (G92) 

The bond is invariant under time-reversal operation, 
therefore we have G p u x+y Gtj- = u x+y , which leads to 
{u x + y , r 3 } = 0. Thus we can choose u x + y to be \ T i ±£t 2 . 
And from G Py , we can obtain u x - y = yr-| — £t 2 . 

c. Third neighbor bond 

Since the G x , G y , Gp xy are trivial, we have 

Pi,i+2x ~ Pi,i+2y — ^ 2 x- (G93) 

It is invariant under time-reversal transformation, 
therefore we have { 112 x^ 3 } = 0. It is also invariant un¬ 
der P yi which leads to [W 2 x>ti ] = 0. Therefore we have 
u 2x = C r i- Note that the third neighbor bond is neces¬ 
sary for the IGG to be Z 2 . 

In conclusion, the ansatz are 


= Ui^i + y — 77T2, 

(G94) 

Ui,i-\-x-\-y = X T l “I" 

(G95) 

Ui,i-\-x—y = xn ^7"2, 

(G96) 

2x = ^i,i-\-2y = CTl, 

(G97) 

ah ± 0 . 

(G98) 


After Fourier transformation, the Hamiltonian be¬ 
comes eiix ± 62 ^ 2 , where 

ei = aj + 2</[cos(2 k x ) + cos(2fc y )] ± 4%cos(fcx)cos(/c y ), 

(G99) 

e 2 = 2?7[cos(A:x) + cos(fc y )] — 4£sin(fcx)sin(/c y ). (G100) 

And the dispersion relation is E± = + e|, it is 

easy to see that when aj is large enough, the spin liquid 
is gapped. 

5. Realization of the AtAt1t 0 t£ state 

PSG elements of this state are: 


G x (i) = G y (i) = r 0 ; 

(G101) 

Gp x {i) = (-) i,+i dT i; 

(G102) 

Gp y (i) = (—) ix+iy iTi[ 

(G103) 

Gp xy {i) = t 0 ; 

(G104) 

GrW = ( 

(G105) 
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The analyze of this state is in much the same way as 
that of the ^AtAtAt 0 t 3 case, here we will only show 
the results: 


7. Realization of the Z 2 AT+T+T 0 t°_ state 
PSG elements are: 


= + y = TjT^y 

(G106) 

Ui^i+x+y = X T 1 + £ T 2i 

(G107) 

Ui t i+ x —y = X T 1 — £ t 2 , 

(G108) 

2 x = — C^"l 1 

(G109) 

a o ^ 

(G110) 

The mean-field Hamiltonian is eiTi + e2T2 + £3X3, where 

ei = aj + 2C[cos(2A; x ) + cos(2A:y)] + 4\cos(A: 2 ,)cos(A:y), 


(Gill) 

£2 = -4^sin(A: a ,)sin(A:y), 

(G112) 

e 3 = 2 g[cos(k x ) + cos(k y )\. 

(G113) 

The energy dispersion is E± = ± y/ef + el 

+ e§. When 

Uq is sufficiently large, the dispersion is 
gapped. 

necessarily 

6. Realization of the Z2At+t < It 0 t p 

state 

PSG elements are: 


G x (i) = G y (i) = t°, 

(G114) 

Gp*(i) = G Py {i) = Gp^ = T\ 

(G115) 

G r {i) = it 3 . 

(G116) 

The term OqT; is invariant under G t 1 ~, therefore we 

have {ctgT/, T3} = 0, hence ctg’ 2 ^ 0 and Og = 

0. 

The gauge transformations G x , G y , G Pj; , Gp y , G Pxb are 
all trivial, therefore the only constraint comes from time- 

reversal transformation. 


From Gr(i)uijGr{jy = —Uij, we have{t%-,r 3 } = 0. 
We can therefore write down a symmetry-allowed ansatz 

which realizes a gapped Z2 spin liquid 


Ui^i+x = Ui,i-\-y = TfiTi + TI2T2 5 

(G117) 

Ui,i+x±.y = xn, 

(G118) 

2 x = ^i,i-\- 2 y = 

(G119) 

1,2 n 

a 0 ^ 0. 

(G120) 


The mean-held Hamiltonian is eiTi + e2T2, where 


G x (i) = G y (i) = t°, (G123) 

G Px (*) = G Py (i) = G Pxh (i) = r°, (G124) 

G r (i) = (G125) 

We can prove that this state is necessarily gapless 

which is protected by the PSG symmetrjC^. Since G x 
and G y are trivial, the mean-held ansatz is manifestly 
translational invariant, therefore we can write down the 
mean-held Hamiltonian as 

if(k) =e' i (k)r /i , (G126) 


where /x = 1,2,3. 

The condition G'j-(i)uijGj-(j) = — u^j , when trans¬ 
lated into k space, becomes g'j-H(k x ,k y )gj- = —H(k x + 
7r, fcy + 7r), gj- = To- Therefore we have 

k y ) = -e /i (fc x . + 7r, k y + 7r). (G127) 


And from Gp^H{k x , k y )G^ P = H(k y , k x ) we have 


£^{k X i ky) £^ (A:?/, k x ), 


(G128) 


Eq. (G127) and Eq. ( G128| ) indicate that the Hamil¬ 
tonian is gapless along the line (k x , k x + tt). On the one 
han d, we have e^{k x ,k x + 7r) = — e^(k x + 7r, k x ) from 
Eq. (G127). On the other hand, we have e fl (k x , k x + w) = 


(k x +ir,k x ) from Eq. (G128I. Therefore e^(k x , k x + 7r) 


must vanish for any k x . indicating that the Hamiltonian 
is gapless along this line. Similarly, we can prove that the 
Hamiltonian is gapless along the line (k x , 7r — k x ). Thus, 
we have completed our proof that the Z 2 At p t p t q t® state 
is necessarily gapless. 


8. Realization of the ZzAt+t+t 0 T- state 
PSG elements are: 

G x (i) = G v {i) = G Px (*) = G Py (i ) = Gp is (*) =(E?129) 
G r (i) = (_)<.+‘»* T 3. (G130) 

The term a l 0 Ti is invariant under GpT, therefore we 
have {uqT/, 73} = 0, hence a Q ’ ^ 0 and clq = 0. 


We can label site i with its parity (—1 Then 

£1 = °o + 2771 [cos(A: a ;) + cos(A; 1/ )] + 4ycos(A: a ;)cos(A;|^1121)f or a k onc i Ui ^ connecting even site i with odd site j, we 
£2 = flo + 2?72[cos(A: a ;) + cos(A;.y)] + 2£[cos(2A; x ) + cos(0£y^2)have [ Uij,r 3 ] = 0. And for a bond Uij connecting two 

even sites or two odd sites, we have {tty, r 3 } = 0. Since 
The dispersion relation is E± = ± \/e 3 + e|. When u ij becomes u p under Gf = ( P xy P y ) 2 and translation, 
aJ,aQ are sufficiently large, the dispersion is necessarily we have an additional constraint Uij = Uji = it A for any 
gapped. bond Uij. 
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Therefore we could write down a mean-field ansatz 


which realizes a gapped Z 2 spin liquid 


= = TjTS) 

(G131) 

Ui,i+x±y = X'Tb 

(G132) 

2x = = ^2 ? 

(G133) 

1,2 / n 
a 0 0 . 

(G134) 


The mean-field Hamiltonian is therefore e\T\ + t- 2 ji + 


£ 3 X 3 , where 

€1 = a}) + 4xcos(k x )cos(k y ), (G135) 

£2 = ccq + 2£[cos(2 k x ) + cos(2fc y )], (G136) 
e 3 = 2r/[cos(k x ) + cos(ky)}. (G137) 

The energy dispersion is E± = ± V e i + e 2 + e 3 - When 
aj and ag are sufficiently large, the dispersion is neces¬ 
sarily gapped. 
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